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ABSTRACT 


This  thesis  extends  the  concept  of  conditional  feedback  to 
sampled  data  systems. 

Dead-beat  controllers  were  used  to  improve  the  disturbance- 
output  response.  The  dead-beat  controllers  were  designed  by  the 
rep-op  method. 

Linear  and  non-linear  sampled  data  systems  were  studied. 

The  optimisation  of  oscillatory  systems  by  the  use  of  conditional 
feedback  and  the  dead-beat  controllers  was  also  studied. 
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CHAPTER  I 


INTRODUCTION 

A  conditional  feedback  system  is  a  system  in  which  the  feedback  path 
is  effective  only  when  there  is  a  disturbance  at  the  plant.  In  the 
absence  of  a  disturbance,  the  conditional  feedback  system  operates  as  an 
open  loop  system. 

The  concept  of  conditional  feedback  was  originally  put  forth  by 
Lang  and  Ham1.  They  discussed  the  continuous  data  conditional  feed¬ 
back  systems  in  their  paper. 


The  basic  configuration  of  a  continous  data  conditional  feedback 
system  is  shown  in  fig  1.1.  Fig  1.2  shows  the  corresponding  conven¬ 
tional  feedback  system. 


Fig.  1.2  -  CONVENTIONAL  FEEDBACK 
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62(8)  is  the  Laplace  transform  of  the  plant. 

G2(S)  is  the  Laplace  transform  of  a  compensator  to  modify  the 
response  of  the  plant,  if  desired. 

H(S)is  the  Laplace  transform  of  the  sensing  device. 

Gc(S)  is  the  Laplace  transform  of  the  compensator  to  modify  the 
disturbance-output  response.  It  will  be  noted  later  that  the 


input-output  response  is  independent  of  GC(S). 

M(S)  is  the  Laplace  transform  of  the  model  to  be  designed. 


U(S)  is  the  disturbance,  R(S)  the  input  and  C(S)  the  output. 
The  Laplace  transform  of  the  output  can  be  written  as 


C(S)  =  _ U(S) _ 

1  +  G1(S)Gc(S)H(S) 

=  _ U(S) 

1  +  G1(S)GC(S)H(S) 


+  {gp( S)  +  M(S)Gc(S)}  G 1  (s)R(s) 

1  +  G1(S)GC(S)H(S) 

1  Gc(S)) 

+  LI  +  M(S)gp(S)JG1 (S)Gp(S)R(S) 

1  4-  G1(S)GC(S)H(S) 


1.1 


Let  M(S)  =  G1(S)G2(S)H(S) 


C(S)  =  _  U(S) _ 

I  +  G1(S)GC(S)H(S) 


+  {l  +  Gt (s)gc(s)h(s)1gt (S)Gp(S)R(S) 

1  +  G1(S)GC(S)H(S) 


1.2 


-h 


=  _ U(S) _  +  R(S)G1(S)G2(S) 

1  +  G1(S)GC(S)H(S) 

C(S)|  =  G1(S)G2(S) 

R^U(S)  =  0 


1.3 

1.4 


C(S)|  =  _ l _  1.5 

U(S)[R(S)  =  0  1+  G1(S)GC(S)H(S) 

Equation  1.4  shows  that  if  the  model  is  given  by  equation  1.2,  the 


input-output  response  is  independent  of  the  feedback  loop.  The  distur¬ 
bance  -  output  response  is  determined  by  the  feedback  loop  and  hence  it 

can  be  designed  to  be  different  from  the  input-output  response.  This  is 
significant  from  the  design  point  of  view.  In  the  design  of  a  conven¬ 
tional  feedback  system  a  compromise  has  often  to  be  made  between  the 
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input-output  response  and  the  disturbance  -  output  response.  They  are 
basically  of  the  same  nature.  This  problem  is  largely  eliminated  by 
the  use  of  conditional  feedback  as  seen  above. 

There  are,  of  course,  other  methods  of  achieving  the  same  e.g. 
by  the  use  of  auxiliary  feedback  loops  or  the  use  of  pre-filters^^, 
but  the  above  method  is  the  most  straight  forward  and  requires  less 
additional  components. 

Another  advantage  arises  from  the  fact  that  the  system  operates 
as  an  open  loop  as  far  as  the  input  signal  is  concerned.  Thus  non- 
linearities  in  the  plant  can  be  tolerated  without  causing  instability. 
Similar  non-linearities  are  built  in  the  model  for  this  purpose. 

If,  however,  non-linearity  in  the  plant  G^S)  causes  instability  in 
the  feedback  loop  when  there  is  a  disturbance,  then  the  compensator  has 
to  be  so  designed  that  the  system  is  stable  in  the  presence  of  a  dis¬ 
turbance  also.  This  modification  in  the  compensator  will  not  affect 
the  input-output  response.  If  it  is  feasible  to  place  the  non-linearity 
in  the  branch  of  the  transfer  function  G2(S) ,  then  the  non-linearity 
will  be  out  of  the  feedback  loop  and  hence  no  problem  of  instability 
will  arise  due  to  it.  The  input-output  response  will  not  be  affected 
by  placing  the  non-linear  plant  in  the  branch  of  the  transfer 
function  G2(S) . 

Another  advantage  of  operating  as  an  open  loop  is  that  low  power 
gains  are  required. 

The  open  loop  system  is  also  advantageous  when  there  are  trans¬ 
port  delays  in  the  system.  These  time  delays  are  often  a  cause  of 
instability  in  closed  loop  systems. 
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The  model  mentioned  above  whose  transfer  function  is 
G1(S)G2(S)H(S)  will  be  simulated  from  simple  R-C-L  elements  and  will 
not  be  subject  to  changes  except  for  the  tolerance  of  the  components. 

The  purpose  of  this  thesis  is  to  extend  the  idea  of  conditional 
feedback  to  sampled  data  systems  and  explore  the  possible  improve¬ 
ments  by  the  use  of  digital  controllers. 

A  secondary  object  of  the  thesis  is  to  establish  the  effective¬ 
ness  of  the  Repetive  Operation  Method  (known  in  short  as  the  Rep-op 
method,  see  appendix  A),  to  design  dead-beat  controllers  for  linear 
and  non-linear  systems.  This  new,  simple,  direct  and  convenient  method 
is  due  to  Professor  Y.  J.  Kingma,  who  has  also  designed  the  equipment 
that  is  used  for  designing  the  dead-beat  controllers  by  this  method. 

Most  of  the  design  of  digital  controllers  and  the  study  of  the 
response  of  systems  to  input  signals  and  disturbances  was  carried  out 
on  the  analog  computer.  The  entire  work  has  been  divided  into  four 
chapters  of  which  the  present  chapter  deals  with  an  introduction  to 
the  thesis. 

Chapter  II  deals  with  linear  sampled  data  systems.  It  was  observed 
here  that,  for  linear  systems,  it  is  possible  to  simulate  a  continuous 
data  model.  The  continuous  data  model  worked  well  for  systems  having 
no  inter-sample  ripple.  In  a  system  with  inter-sample  ripple,  there 
will  be  an  error  at  the  branch  Gq(S)  during  a  sampling  period.  With 
the  use  of  conditional  feedback,  there  are  possibilities  of  removing 
or  at  least  appreciably  reducing  the  inter-sample  ripple  as  discussed 
in  chapter  IV.  This  is  achieved  by  modifying  the  model. 

Continuous  data  compensators  were  first  tried  and  they  worked  well. 
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The  feedback  loop  in  one  of  these  cases  was  unstable  because  of  high 
gain.  It  was  stabilized  by  designing  the  compensator  GC(S)  by  the  Bode 
Plot  Method  (Bi-linear  transformation  method) . 

It  was  observed  that  a  system  designed  for  a  step  input  worked 
well  for  a  ramp  input  too.  This  was  also  proved  analytically. 

Dead-beat  digital  controllers  were  then  tried.  The  insertion 
of  these  controllers  did  not  affect  the  input-output  response  but  it 
removed  the  effect  of  disturbances  on  the  output  in  an  optimum  manner. 

Chapter  III  deals  with  non-linear  sampled  data  systems.  In  this 
case  a  simple  continuous  model  cannot  be  found  because  there  is  no 
transform  for  the  non-linearity.  The  model  includes  a  sampler  and  a 
zero  order  hold. 

The  dead-beat  controllers  were  designed  by  the  rep-op  method.  This 

chapter  and  the  next  establish  the  convenience  of  the  use  of  rep-op 

method  over  other  methods.  Just  like  the  State  Variable  Method,  the 

rep-op  method  has  the  advantage  that  there  is  no  inter-sampling  ripple, 

but  it  is  much  quicker  than  the  former.  The  transfer  function  of  the 

plant  used  for  a  saturation  type  of  non-linearity  was  1  and  the 

S(S+1) 

non-linearity  was  located  in  the  positions  shown  in  fig.  1.3. 


Fig.  1.3  a 
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Fig.  1.3  b 


The  dead-beat  controller  designed  for  ramp  disturbances  was  modified 
by  trial  and  error  to  give  a  good  response  to  step  disturbances  also. 

The  non-linear  systems  were  studied  for  all  possible  types  of  internal 
and  external  disturbances.  The  response  was  good  in  most  of  the  cases. 
The  effect  of  change  of  saturation  limit  of  the  non-linearity  on  the 
output  was  also  studied. 

The  system  used  for  hysteresis  type  of  non-linearity  was 


(SH) 

The  simulation  of  hysteresis  non-linearity  on  the  analog  computer  was 

carried  out  as  shown  in  appendix  B.  The  hysteresis  non-linearity  was 

placed  in  the  branch  of  the  sensing  device. 

Chapter  IV  deals  with  a  second  order  oscillatory  system  with  a 

2 

transfer  function  of  the  form  ^  f  ^  \  .  Normally  a 

S2-  -f  2  5  U3v\  -V  tOiy^ 

system  like  this  will  have  overshoot  and  oscillations  before  settling 
down  to  its  steady  state  value,  but  it  was  found  that  if  a  proper  damping 
ratio  is  chosen  for  the  model,  an  optimum  response  can  be  obtained. 

It  was  observed  here  that  the  smaller  the  sampling  period  of  the  con¬ 
troller  the  better  the  response. 

The  optimum  value  of  the  damping  ratio  for  the  model  was  chosen 
by  trial  on  the  analog  computer.  The  overshoot  of  an  oscillatory  sys¬ 
tem  having  a  damping  ratio  of  0.5  could  be  reduced  to  almost  zero 
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percent  with  T  =  0.2  sec.  For  a  damping  ratio  of  0.3,  the  overshoot 
could  be  reduced  to  3%  with  T  =  0.2  sec. 

A  conventional  feedback  system  having  a  dead-beat  controller  in 
it  is  very  sensitive  to  parameter  changes  in  the  plant.  The  con¬ 
ditional  feedback  system,  optimised  for  the  output,  and  the  cor¬ 
responding  conventional  feedback  system  with  dead-beat  controller 
were  compared  for  the  same  changes  in  parameters.  It  was  found  that 
the  performance  of  the  conditional  feedback  system  was  far  superior 
to  that  of  the  conventional  feedback  system. 
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CHAPTER  II 


This  chapter  deals  with  linear  sampled  data  systems. 
2.1  TYPE  ZERO  SYSTEM 

Ki 


A  plant  with  the  transfer  function 


was  chosen. 


(S  +  1) 


Its  response  was  modified  by  the  transfer  function 
in  fig.  2.1. 


Ks 


as  shown 


(S  +  2) 


0.861  (  -  S  +  1. 16)^2 
(S2  +  3S  +  2) 


Ki  =  1 
K2  -  1 


Fig.  2.1 

H(S)  was  assumed  to  be  unity.  The  model  can  be  derived  as  follows: 

Ki 
*  s+T 

C(S)  =Gho(S)G1(S){R(S)G2(S)]*t 
C*(S)  =[Gho(S)G1(S)}  *[r(S)G2(S)}  * 

C(Z)  =  GhoG^Z)  x  RG2(Z) 

T  =  1  sec. 


^  *  Shows  the  sampled  signals 
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Now  Gho  (8)02(8) 


,  -  TS  Kn 

1  -  e 

S  X  (S+l) 


and  GhoG1(Z) 


0 . 632K-L 
(Z  -  0 „ 368) 


1  K2 

R(S)xG2(S)  =  -  x  s  (S+2) 


0 . 432K2Z 

RG2(Z)  =  (Z-l)  (Z-0. 136) 


0 . 432K2Z  0.632K! 

C(Z)  =  (Z-l)  (Z-0.  136)  X  (Z-0 . 368) 

0.2725K2K2Z 

=  (Z-l)  (Z-0 . 136)  (Z-0„ 368) 


Splitting  ^ 


into  partial  fractions 


z 


0.2725^2 


1.833  4. 98 

(Z-l)  (Z-0. 136) 


6,82  1 

(Z-0. 368) J 


Taking  the  inverse  2-  transform: 


C(S)  =  0 . 2725K-lK2 


f  1-833 

L  S 


+ 


4,98 

S+2 


6.82 

S+l 


0. 861K1K2(-S+1. 16) 
S  (S+l)  (S+2) 


This  must  also  be  the  Laplace  transform  of  the  output  of  the  model 

0 . 861K1K2(-S+1 . 16) 

R  (S)  xM(S)  = 

S(S2+3S+2) 


R(S)  =  i 

0.861K1K2(-S+1. 16) 
M(S)  =  . . .  "  . . 


(S2+3S+2) 


-X 


10 


The  signal  flow  graph  for  the  model  is  shown  in  fig.  2.2. 


Fig.  2.2 


The  analog  computer  diagram  of  the  conditional  feedback  arrangement  is 
shown  in  fig.  2.3  ref. 


Fig.  2.3 

The  response  of  the  system  to  a  step  input  is  shown  in  fig.  2.4. 
This  is  with  GC(S)  as  1.  GC(S)  was  then  assumed  to  be  K  and  Routh- 
Hurwitz  criterion  applied  to  study  the  stability  of  the  feedback  loop. 
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Fig.  2.4. 


The  characteristic  equation  is 


1  +  KGhoGx (Z)  =  0 


!  +-°-632K  ■  =  0 

1  +  (Z-0.368) 


Z  -  0.368  +  0.632K  =  0 


Applying  the  bi-linear  transformation 


r+1  .  i j 

Z  =  ““7  yields, 
r-1 


r+1 

r-1 


-  0.368  +  0.632K  =  0 


r  = 


0.632K  -  1.368 

0.632K  +  0,632 


The  root  r  must  lie  in  the  left  half  r-plane  for  the  system  to  be  stable 


Hence,  the  system  is  unstable  for 
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v  1.368 
>  0.632 


K  >2.17 


This  value  was  verified  on  the  analog  computer  as  2.15. 

The  feedback  loop  for  the  disturbance  is  shown  in  fig.  2.5. 


Fig.  2.5 

The  steady  state  value  of  C(S)  due  to  a  step  disturbance  U(S)  can 
be  calculated  as  follows; 

C(S)  =  E  (S) 

E(S)  =  U(S)  -  KGho(S)  G1(S)C*(S) 

U*(S) 


C*(S)  = 


1  +  K{Gho(S)G1(S)j  * 


C(Z)  = 


M 


1  +  KGhoGj.CZ) 
Applying  the  final  value  theorem 


“1 

Lim.  C*(t)  =  Lim.  (1  -  Z  )  U(Z) _ 

t->~  Z->1  1  +  KGhoGi  (Z) 


=  Lim. 

Zo  1 


1  +  K 


0.632 

(Z  -  0.368) 


1 

1  +  K 
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This  shows  that  the  larger  the  value  of  K,  the  smaller  will  be  the 

steady  state  effect  on  the  output  due  to  a  step  disturbance. 

K  was  taken  as  5  and  then  a  digital  compensator  was  designed  to  make 

(2) 

the  system  stable.  The  bi-linear  transformation  method  was  used 
as  follows: 


GhoG1(Z) 


0.632K 

(Z  -  0.368) 


Applying  W  transformation 


1  +  W 

1  -  W 


GhoG^W) 


0.632K 
1  +  W 

1  -  W 


0.368 


0.632(1  -  W)K 

0.632  +  1.368W 


K  (1  -  W) 

1  +  2.17W 


GhoG^jw  ) 
w 


K(1  -  J»w) 

1  +  2. 17 jw 

w 


The  Bode  diagram  for  1(1  -U>j  w)  plotted  in  fig.  2.6.  Phase 

(1  +  2.17  jww) 

lag  compensation  was  chosen.  Phase  lead  compensation  increases  band¬ 
width  and  is  also  not  so  effective  if  the  phase  lag  increases  rapidly  in 
the  vicinity  of  the  gain  cross  over.  The  increase  in  the  bandwidth  is 
not  desirable  in  sampled  data  systems.  The  sampling  switch  produces 
harmonics  at  its  output.  If  the  bandwidth  is  increased,  the  harmonics 
will  be  less  suppressed  and  the  response  will  not  be  so  smooth. 

The  new  gain  cross-over  was  chosen  as  0.6.  The  uncompensated  system 
has  a  gain  of  11.5  db  at  this  frequency. 
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db 


The  phase  lag  compensator  is  of  the  form 


Gc(W) 


1  +  aTW 
1  +  W 


Then  20  log  a  =-11.5 
,  1  =  0.575 

log  7 

a  =■  0.266 

Placing  the  upper  corner  frequency  one  decade  below  the  new  cross¬ 
over  frequency. 


1  =0.06 
a.'t 


^  =  62.65 
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GC(W) 


1  +  16, 7W 

1  +  62.65W 


putting  W  = 


Z  -  1 
Z  +  1 


yields , 


GC(Z) 


0.278 


1  -  0.887Z 

“  1 


1  -  0.968Z" 

Gain  margin  of  this  system  was  4~  db  and  the 
Response  to  a  step  disturbance  with  this  controller 
2.7  -  1*.  The  settling  time  was  observed  to  be  too 
compensator . 

The  new  gain  cross-over  was  then  chosen  as  0.8 


phase  margin  100°. 
is  shown  in  fig. 
long  with  this 

and  D(Z)  was 


D(Z)  =  0.333 


"l 

1  -  0.8525Z 

"l 


1  -  0.952Z 
Gain  margin  =  2  db 

Phase  margin  =  81° 

The  response  to  a  step  disturbance  is  shown  in  fig.  2.7  -  2. 

Next  K  was  chosen  as  10  and  with  the  new  gain  cross-over  as  0.8. 


D  (Z)  =  0.169 


"1 

1  -  0.8525Z 

“l 


1  -  0.975Z 
Gain  margin  =  2  db 

Phase  margin  =  816 

A  free  integrator  was  then  tried  as  the  compensator.  Analytically 
the  effect  of  a  step  disturbance  on  the  output  with  this  compensator 
should  be  zero  in  the  steady  state.  The  system  was  tried  on  the  analog 
computer  and  it  showed  a  much  too  oscillatory  response  to  a  step  dis¬ 
turbance  . 


*  dash  shows  the  curve  number  on  that  figure. 


' 
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Fig.  2.8 
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A  dead-beat  controller  was  then  tried.  This  is  a  type  zero  sys¬ 
tem  and  therefore  an  integrator  has  to  be  used  with  the  controller,  so 
that  when  the  input  to  the  controller  goes  to  zero  the  integrator  should 
hold  the  output  value (see  appendix  A) . 

The  feedback  loop  in  this  case  is  shown  in  fig.  2.9. 


a  simple  case  like  this,  the  dead-beat  controller  can  be  derived  step 

K 


by  step  as  shown  below.  A  general  transfer  function 


(S  4-  C) 


is 


considered. 

Sampling  period  =  T  sec. 

The  signals  at  different  points  are  shown  in  fig.  2.10. 

The  digital  controller  can  be  represented  as  shown  in  fig.  2.11 
The  disturbance  is  a  unit  step.  Therefore,  the  input  to  the  con¬ 
troller  at  the  first  sampling  instant  is  1.  This  multiplied  by  the  gain 
'a'  appears  at  the  input  to  the  integrator.  Nothing  will  be  applied  at 
the  input  to  the  plant  at  the  first  sampling  instant.  During  the  first 
sampling  period,  the  output  at  the  integrator  will  rise  to  'aT*  and  at 
the  next  sampling  instant  a  step  input  of  ’aT'  will  be  applied  to  the 

plant.  _Ct 

.  ,  aT  K  (1  -  e  ) 

The  response  of  the  plant  is  given  by  . — 

C 


. 
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0  T  2T  3T 
(a)  Disturbance  input 


0  T  2T  3T 

(b)  Output  of  the  plant 


(d)  Output  of  the  integrater 


(e)  Input  to  the  plant 
Fig.  2.10 

For  dead-beat  response,  this  should  be  unity  at  the  end  of  one 
sampling  period 

-  CT 

aT  K  (1  -  e  ) 

C 


1 


< 
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O 


Fig.  2.11 


C 

a  =  — - — — — — — 

-  CT 

T  K  (1  -  e  ) 

Next  the  value  of  the  gain  ' b 1  is  to  be  found.  At  the  sampling 
instant  t  =  T  the  negative  feedback  of  the  time  delay  will  cancel  the 
unit  disturbance  input  to  the  controller.  A  signal  has  now  to  be  applied 
to  the  input  of  the  integrator,  so  that  at  the  sampling  instant 
t  =  2T,  the  output  of  the  integrator  is  unity.  The  value  of  ' b  *  will  be 


(1  -  a) 

T 

f  C  1 

°r  b  =  1  -  - ITT  T 

1  TK(1  -  e  )J 

After  that  the  output  of  the  integrator  will  stay  at  unity  and  the  input 
to  the  controller  at  zero. 

In  this  particular  case 


K=l,  T  =  l,  C  =  1 
a  =  1.582,  b  =-0.582 


D  (Z) 


"1 

1.582  -  0.582Z 

“l 

1  +  Z 


The  response  of  the  system  to  a  step  input  with  this  dead-beat 
controller  as  the  compensator  is  shown  in  fig.  2.8-1.  A  step  distur- 
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bance  was  applied  at  the  output  of  the  plant  when  the  response  to  the  step 
input  had  settled  down  to  its  steady  state  value.  The  output  of  the 
plant  was  brought  back  to  its  original  value  in  two  sampling  periods. 

The  process  of  removal  of  the  disturbance  will,  start  from  the  sam¬ 
pling  instant  immediately  following  the  application  of  the  disturbance. 
Thus  the  actual  time  of  removal  of  the  disturbance  can  be  anywhere 
between  two  and  three  sampling  periods.  This  additional  time  Ta  is 
given  by  the  following  inequality 

0  <  Ta  <  T 

This  comment  about  the  additional  time  will  hold  for  all  other 
cases  to  be  studied. 

The  step  disturbance  was  then  applied  at  the  point  ' 1'  in  fig.  2.1 
and  the  response  is  shown  in  fig.  2.8  -  2.  This  disturbance  may  also  be 
considered  as  a  variation  in  the  parameter  or  the  gain  of  the  transfer 
function  G2(S)  ,  because  that  will  have  a  similar  effect  at  the  point  '1'. 

The  response  to  a  step  disturbance  at  'm'  in  fig.  2.1  is  shown  in 
fig,  2.8  -  3.  It  may  be  noted  here  that  a  disturbance  at  'm'  gives  rise 
to  an  equal  disturbance  at  the  output,  hence  care  will  have  to  be  taken 
that  the  transfer  function  H(S)  is  not  subject  to  large  disturbances. 

An  analog  integrator  in  the  controller  is  undesirable,  because  if 
there  is  any  drift,  it  will  start  building  up  at  the  integrator.  More¬ 
over,  the  use  of  the  analog  integrator  requires  an  additional  sampler 
and  a  zero  order  hold.  The  zero  order  hold  and  the  analog  integrator  may 

be  combined  with  the  controller.  The  Z-transform  of  a  zero  order  hold 

"l 

and  an  integrator  is  Z  .  The  digital  integrator  is  shown  in 

(i  -  z1) 


fig.  2.12. 


V 
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Fig.  2.12  Digital  integrater 


Modified  D(Z)  =  X,1-582  ~  8-582Z  1 

(1  +  Z  ) 


x 


(1  -  Z  ) 


"1  “2 

D(Z)  =  1.582Z  -  0.582Z 

m  (i  -  z2) 


This  was  tried  as  shown  in  fig.  2.13  -  1.  It  worked  well. 


Internal  variations  in  the  plant  such  as  gain  variations  and  time 


K-, 


constant  variations  were  studied.  For  the  transfer  function 


(S  +  1)’ 


an  instantaneous  gain  variation  or  a  time  constant  variation,  will  be 
equivalent  to  the  sum  of  step  and  exponential  disturbances  at  the  out¬ 
put  of  the  plant.  The  dead-beat  controller  designed  for  step  distur¬ 
bances  was  found  to  behave  well  for  these  combined  disturbances  also. 
The  results  are  shown  in  table  2.1 

These  changes  in  the  parameter  were  affected  instantaneously  by 
switches.  In  actual  practice  the  change  is  likely  to  occur  gradually 
and  therefore  the  output  will  not  be  disturbed  too  much. 

The  unevenness  in  the  transient  response  is  due  to  the  fact  the 
input  to  the  plant  and  the  model  is  not  applied  at  the  same  moment. 
Signal  to  the  plant  is  applied  at  the  sampling  instant  that  follows  the 


'  f  - 
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Fig.  2.14 


\ 
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application  of  the  input.  The  model  does  not  contain  a  sampler,  so 
the  signal  is  applied  to  it  immediately  after  the  application  of  the 
input . 


TABLE  2.1 


Fig. 

Parameter 

Original 

New 

Changed 

Value 

Value 

2.13  -  2 

k2 

1 

1.1 

2.13  -  3 

k2 

1 

1.5 

2.14  -  1 

Time  constant 

1 

1 

(S  +  1) 

(S 

+  1.1) 

2.14  -  2 

Time  constant 

1 

1 

(S  +  1) 

(S 

+  1.5) 

Next  a  ramp 

input  was  applied 

to  the  above  system. 

It 

removed  the 

external  disturbances  as  before. 

This  is  shown  in  fig.  2 

.15 

».  The 

application  of  a 

ramp  input  can  be 

analyzed  analytically 

as 

shown  below 

The  response 

of  the  model  to 

a  ramp  input  will  be 

V(S)  =  R(S)  x  M(S) 

0 . 861K-L  (-S  +  1.16)  __  _  . 

_ _ _  5  ^2.  ” 

S2(S  +  1)  (S  +  2) 

Splitting  into  partial  fractions  yields, 


V(S) 

Taking  the 


v  / J_  0.25  1.86 

1(v2$2  S  S  +  1 

inverse  Laplace  transform 


0.68  \ 
S  +  2/ 


V(t)  =  0.5  Kxt  -  0.25  Kx  U(t)  +  l^K^  -  O^S^e^ 

The  response  of  the  plant  at  the  sampling  instants  can  be  derived 


as  follows: 
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G2(s)  G1(S) 


Fig  2.15  (a) 


The  open  loop  system  is  shown  in  fig.  2.15  (a). 


[  i  I  *  -TS  Ki 

C(S)  -  \  R(S)  x  c  ^  0  }  x  (1  -  e  )  x 


S  +  2 


S(S  +  1) 


j  1  {*  j  -TS  Ki  1  * 

C*(s)  =f  R(S)  X  X  f(l  -  e  ) 


S  +  2 


S(S  +  1) 


c (Z)  =  RG2(Z)  x  (1  -  z  )  ?J  Kl 

0  ls(s  +  i)J 

RG  (Z)  =  0 . 5Z (0 , 568Z  +  0.296) 

(Z  -  1)2(Z  -  0.136) 

0 . 5Z (0 . 568Z  +  0.296)  0.632% 

C(Z)  =  -  x  - 


(Z  -  1)2(Z  -  0.136) 


(Z  -  0.368) 


Dividing  C (Z)  into  partial  fractions 
Z 


Z 


0.316% 


1.585 

(Z  -  1)  2 


3.3  5.45 

Z  -  1  Z  -  0.368 


2.155 

Z  -  0.136 


Taking  the  inverse  Z-transform 


- 1  -2t 

C*(t)  =  0.5% t  -  1 . 04K1U( t)  +  1.72Kie  -  0.68Kie 


Comparing  the  expressions  for  C*(t)  and  C(t),  it  is  found  that  the 
two  will  differ  by  a  step  signal  in  the  steady  state.  This  comparison 
is  at  the  sampling  instants  only,  but  if  there  is  no  inter-sample  ripple, 
the  comparison  will  hold  at  all  the  points  of  the  output.  This  step 
signal  will  be  treated  as  a  disturbance  in  the  feedback  loop  and  hence 


. 

> 
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Fig 


Fig.  2.16 
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compensated  for. 

The  sampler  was  then  placed  outside  the  feedback  loop  i.e.  in 
the  branch  of  the  transfer  function  G2(S) .  This  will  be  useful  when 
it  is  desired  that  the  sampling  period  should  not  have  any  effect  on 
the  stability  of  the  feedback  loop  and  also  when  it  is  not  desired 
to  use  expensive  digital  controllers.  The  placing  of  the  sampler  in 
the  branch  of  G2(S)  will,  however,  not  always  be  possible.  The 
outputs  were  as  shown  in  table  2.2. 

TABLE  2.2 


Fig. 

Remarks 

2.16  -  1 

GC(S)  = 

1 

S 

2.16  -  2 

GC(S)  = 

10 

2.16  -  3 

Gc (S)  = 

100 

2.2  TYPE  ONE  SYSTEM 

The  system  studied  in  this  case  is  shown  in  fig.  2.17. 


S(S  +  1)  (S  +  2) 


Fig.  2.17 


, 
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The  model  can  be  derived  as  follows 
C  (S)  =  [r(S)G2(S)|  *  [Gho(S)G1(S)j 
C*(S)  =[r(S)G2(S)]  *  [Gho(S)G1(S)j  * 
C(Z)  =  RG2(Z)  x  GhoGi (Z) 


-TS 

1  -  e  K-l 


Gho(S)G1(S) 

s  A  S(S  +  1) 

GhoG1(Z)  = 

0.368K1(Z  +  0.713) 

(Z  -  1)  (Z  -  0.368) 

R(S)G2(S)  = 

k2 

S(S  +  2) 

RG2(Z)  =  — 

0 . 432K2Z 

(Z  -  1)  (Z  -  0. 136) 


0 . 159K7 K2  Z  (Z  +  0.713) 

C(Z)  =  - 

(Z  -  1)2(Z  -  0.136)  (Z  -  0.368  ) 

Dividing  C (Z)  into  partial  fractions 
Z 


3.14 


6.77 


4.9 


c  (z)  =  kxk2 


ii.: 


+ 


■(Z  -  l)2  (Z  -  1)  (Z  -  0.136)  (Z  -  ( 


Taking  the  inverse  Z~  transform 


3.14 


6.77 


4.9 


11.7 


C(S)  = 


kxk2 


=  0.2113^2 


S  S  +  2  S  +  1 

(S2  -  3.09S  +  4.77) 

S2(S  +  1)  (S  +  2) 


.368)  J 


This  must  also  be  the  Laplace  transform  of  the  output  of  the  model 


-  - 


■3 
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|  x  M(S)  =  0.211310^2 


(S2  -  3.09S  +  4. 77) 

S(S  +  1)  (S  +  2) 


;,M(S)  =  0. 2113K1K2 


(S2  -  3.09S  +  4.  77) 

S(S  +  1)  (S  +  2) 


The  signal  flow  graph  of  the  model  is  shown  in  fig.  2.18  (a). 

The  computer  diagram  is  shown  in  fig.  2.18  (b) . 

Since  it  is  a  type  one  system  and  it  is  operating  in  the  open  loop 
configuration,  there  will  be  a  ramp  output  for  a  step  input.  This 
suggests  that  a  parameter  variation  in  the  plant  will  be  equivalent  to 
the  sum  of  step,  exponential  and  ramp  disturbances  at  the  output  in  the 
transient  stage.  In  the  steady  state  this  will  be  equivalent  to  the  sum 
of  step  and  ramp  disturbances.  Exponential  disturbances  will  die  out. 

In  this  case  a  suitable  digital  controller  will  be  the  one  which  gives 
a  good  response  to  both,  the  step  and  the  ramp  disturbances. 

A  dead-beat  controller  for  a  step  disturbance  will  give  a  constant 
error  to  a  ramp  disturbance  whereas  a  dead-beat  controller  for  a  ramp 
disturbance  will  give  a  high  overshoot  to  a  step  disturbance.  A  com¬ 
promise  is,  therefore,  to  be  made  between  the  two. 

(2) 

The  combined  synthesis  and  minimal  response  method  was  used  to 
find  the  controller.  This  method  was  found  to  be  quite  convenient  for 
a  case  like  this.  This  method  allows  one  to  design  the  digital  controller 
according  to  the  desired  maximum  overshoot  and  the  settling  time. 

The  maximum  overshoot  and  the  settling  time  are  related  to  the  damping 
ratio  and  the  damping  constant  respectively.  The  requirements  of  the 
proper  damping  ratio  and  the  damping  constant  can  be  met  because  one 
is  free  to  choose  the  roots  of  the  characteristic  equation. 


. 


1 


0.2113K,K 
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The  feedback  loop  can  be  viewed  as  shown  in  fig.  2.19. 


Fig.  2.19 


T  =  1  sec. 


Gp(Z) 


-TS 


1  -  e 


L  S2(S  +  1)  J 

"i  “i 

0.368  Z  (1  +  0.713Z  ) 
(1  -  Z1) (1  -  0.368Z1) 


Since  Gp(Z)  has  one  more  pole  than  it  has  zeros,  to  insure  that 

D(Z)  is  physically  realizable  the  overall  function  M(Z)  /  M(Z)  =  C (Z) 

V  R(Z) 


must  have  at  least  one  more  pole  than  it  has  zeros.  To  satisfy  the 

zero  steady  state  error  requirement  of  the  ramp  response,  (1  -  M(Z)) 

~1  2 

should  contain  the  factor  (1  -  Z  )  .  This  also  satisfies  the  require¬ 
ment  that  (1  -  M(Z) )  must  include  the  poles  of  Gp(Z)  which  lie  on  or 
outside  the  unit  circle  in  the  Z-plane. 


Now  M(Z) 


C(Z)  D(Z)Gp(Z) 
R(Z)  1  +  D(Z)Gp(Z) 


The  poles  and  zeros  of  Gp(Z)  inside  the  unit  circle  are  cancelled 
by  the  poles  and  zeros  of  D(Z) .  Therefore  D(Z)  is  of  the  form 


' 
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D(Z)  = 


M(Z) 


1  -  M(Z) 


x 


Gp(Z) 


=  k. 


(i  -  o^bsz1^  i  +  ai z ^(l  +  a^z1)  - 

(l  +  oxioz^a  -  z^-xi  +  b^z1) - 


For  D(Z)  to  be  physically  realizable,  the  number  of  poles  should  be 
at  least  equal  to  the  number  of  zeros.  In  practice,  the  minimum 
possible  number  of  poles  and  zeros  will  be  chosen  so  that  the  simu¬ 
lation  of  the  digital  controller  is  the  simplest. 

The  open  loop  transfer  function  of  the  compensated  system  is 

O.abSkdZ1  (1  +  a1Z1)(l  +  a2Z1)  - 

D(Z)Gp(2)  =  .  . .  ■■ 

(1  -  Z1)2(l  +  bjZ1) (1  +  b2Zx)  - 


The  lowest  power  in  Z  of  the  system  is  two, 
simplest  form  of  D(Z)Gp(Z)  is 

O.SbSkdZ1  (1  +  a^Z1) 


D(Z)  Gp(Z)  = 


(1  -  z1)2 


In  other  words,  the 


which  has  two  unknowns  kd  and  a±.  The  characteristic  equation  of  the 

0.368kdZx  (1  +  a^Z1) 

system  is  1  H - -  =  0 

(1  -  Z1)2 

1  +  (0.368  kd  -  2X1  +  (1  +  0. 368kd  a2)  Z2  =  0  2.1 


The  damping  ratio  curves  shown  in  fig.  2.20  can  now  be  used  and 
suitable  poles  chosen  for  the  system.  The  poles  were  first  chosen  as 
0.27  +  j  0.32,  to  give  a  damping  ratio  of  0.707.  The  characteristic 
equation  with  these  roots  is 


(Z  -  0.27  -  j  0.54)  (Z  -  0.27  +  j  0.54)  =  0 
1  -  0.54Z1  ^  0.1754Z2  =  0 


2.2 


(.  'j; 

CJ  '  .  •  -0  (  T  1 

.. 


- 


■:L .  r:-  :.c  -  •) 
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Equating  the  coefficients  of  equal  powers  in  equations  2.1  and  2.2. 
0.368kd  -  2  =  -0.54 
kd  =  3  .967 

and  1  +  0.368kda;L  =  0.1754 
ax  =  -0.563 


"l 

(1  -  0.368Z  ) (1  -  0.563Z1) 

D(Z)  =  3.967  - - - r- - 

(1  +  0.713Z  ) (1  -  Z  ) 


“1  "2 

(1  -  0.931Z  +  0.2072Z  ) 

-3.%7  - - - — 

(1  -  0.267Z  -  0.713Z  ) 

Then  the  roots  were  chosen  as  0.5  +  j  0.15  to  give  a  damping  ratio  of 
0.88  and  the  D(Z)  found  was 


■ 

;r  9  '  •  ■  ;  -  <:  .  •*  :  IV  -oJrj  ■;  -V  f  i>  ft  “ 


fii  ■ 
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D (Z)  =  2.72 


"1  “2 

(  1  -  1.0955Z  +  0.268Z  ) 

(1  -  0.287Z1  -  0.713Z2) 


A  dead-beat  controller  for  ramp  disturbances  was  then  designed. 

The  performance  of  these  controllers  was  then  compared.  The  design 

(2) 

procedure  for  the  dead-beat  controller  is  shown  below. 

Since  there  is  a  ramp  input  and  Gp(Z)  has  no  zeros  outside  the  unit 
circle,  M(Z)  will  be  of  the  form 

"1  "2 

M(Z)  =  2Z  -  Z 


D  (Z) 


M(Z)  __1 _ 

1  -  M(Z)  X  Gp (Z) 


Substituting  the  values  of  M(Z) 


and  Gp(Z) 


“1  "2  "1  "1 

(2Z  -  Z  ) (1  -  Z  )  (1  -  0.368Z  ) 

D  (Z)  =  - : - 1 - : - rrj— 

(1  -  2Z1  +  Z2) (0.368Z1) (1  +  0.713Z  ) 

"1  "2 

(1  -  0.868Z  +  0. 184Z  ) 

=  5.44  - - - =S - 

(1  -  0.287Z  -  0.713Z  ) 


These  controllers  were  then  used  in  the  system  of  fig.  2.18  and 
the  results  are  shown  in  tables  2.3  to  2.5. 

1 .  Dead-beat  controller  for  ramp  disturbances 


"1  "2 

(1  -  0.868Z  +  0.184Z  ) 

D (Z)  =  5.44  - — - rj- 

(1  -  0.287Z  -  0.713Z  ) 

T  =  1  sec. 


Max.  overshoot  to  a  step  disturbance  =  1507. 
Settling  time  for  a  ramp  disturbance  =  3  T 


Settling  time  for  a  step  disturbance  =  13  T 


- 


; 
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Fig.  2.22 
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Fig.  2.23 


Fig.  2.24 
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TABLE  2.3 

STEP  INPUT  TO  THE  SYSTEM  =  2  V 


Fig. 

2.21  -  1 
2.21  -  2 


Remarks 

Response  without  disturbance 
the  plant  is  changed  to 
1 

S(S  +  1.5) 


2.21  -  3  A  step  disturbance  of  2  V  is 

applied  at  the  output 

2 .  The  closed  loop  poles  of  the  compensated  system  at 


(0.27  +  j  0.54) 

"1  “2 
(1  -  0.931Z  +  0.207Z  ) 

D (Z)  =  3.967  - — - - — 

(1  -  0.267Z  -  0.71.3Z  ) 

T  =  1  sec. 

Max  overshoot  to  a  step  disturbance  =  1007o 
Settling  time  for  a  step  disturbance  =  8T 
Settling  time  for  a  ramp  disturbance  =  4T 
Step  input  to  the  system  =  2  V. 


Fig. 


2.22  -  1 
2.22  -  2 
2.22  -  3 


TABLE  2.4 

Parameter  Original  New 

Changed  Value  Value 

Response  without  disturbance 
K2  1  1.1 

K2  1  1.3 

K2 


2.22  -  4 


1 


1.5 


! 
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Fig. 

Parameter 

Original  New 

Changed 

Value  Value 

2.22  -  5 

Step  disturbance  of  2  V  at  the  output 

2.23  -  1  &  2 

Time  constant 

1  1 

S(S  +  1)  S(S  +  1.  1) 

2.23  -  3 

Time  constant 

1  1 

S(S  +  1)  S(S  +  1.3) 

2.23  -  4 

Time  constant 

1  1 

S(S  +  1)  S(S  4-  1.5) 

2.23  -  5 

A  step  disturbance 

of  2  V  is  applied  at  the  point 

' 1 '  in  fig.  2.17. 

2.23  -  6 

A  step  disturbance 

of  2  V  is  applied  at  the  point 

'm'  in  fig.  2.17. 

3 .  The  closed  loop  poles  of  the  compensated  system  at 

(0.5  ±  j  0.15) 

“1  “2 
(  1  -  1.095Z  +  0.268Z  ) 

D(Z)  =  2.72  - — - -=y— - — 

(  1  -  0.287Z  -  0.713Z2) 

T  =  1  sec. 

Max.  overshoot  to  a  step  disturbance  =  38% 

Settling  time  for  a  step  disturbance  =  13  T 

Settling  time  for  a  ramp  disturbance  =  8  T 

Step  input  to  the  system  =  2  V 

TABLE  2.5 

Fig.  Parameter  Original  New 

Changed  Value  Value 


2.24  -  1 


Response  without  disturbance 


. 
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Fig 

Parameter 

Original 

New 

Changed 

Value 

Value 

2.24  - 

2 

k2 

1 

1.5 

2.24  - 

3 

Time  constant 

1 

1 

S(S  +  1) 

S(S  +  1.5) 

2.24  - 

4 

Time  constant 

1 

1 

S(S  +  1) 

S(S  +  1.2) 

2.24  - 

5 

A  step  disturbance  of  2  V  is 

applied  at  the 

output . 

Here  again  the  change  of  parameters  was  affected  instan¬ 
taneously  by  switches.  In  actual  practice,  they  are  likely  to 
occur  gradually  and  thus  the  output  will  be  less  affected. 

A  ramp  input  was  then  applied  to  the  system.  It  gave  a  good 
parabolic  response.  The  error  in  the  branch  of  GC.(S)  diminished 
to  zero.  The  controller  effectively  removed  step  disturbances 
at  the  output.  Any  parameter  change  was  equivalent  to  a  para¬ 
bolic  disturbance  and  since  the  controller  had  been  designed  for 
ramp  disturbances,  it  gave  a  constant  error  to  acceleration  dis¬ 
turbances.  If,  however,  ramp  inputs  are  generally  expected  into 
the  system,  the  controller  could  be  designed  for  parabolic  dis¬ 
turbances  and  then  modified  to  give  a  good  response  to  ramp  and 


step  disturbances. 


. 


CHAPTER  III 


This  chapter  deals  with  non-linear  sampled  data  systems. 

3.1  SATURATION  TYPE  NON-LINEARITY  BEFORE  THE  PLANT 

The  transfer  function  of  the  plant  chosen  in  this  case  was 

_ 1  .  The  conditional  feedback  arrangement  is  shown  in  fig. 

S(S  +  1) 

3.1. 


T  =  1  sec. 

Fig.  3.1 


Step  input  into  the  system  was  considered.  In  this  case  the 
model  was  simulated  exactly  similar  to  the  actual  system  i.e.  it 
included  a  sampler,  a  zero  order  hold  and  the  non-linearity.  A 
continuous  model  could  not  be  found,  because  there  is  no  transform 
for  the  non-linearity. 

The  dead-beat  digital  controller  was  designed  by  the  rep-op 
method  for  a  ramp  disturbance  of  2  V/sec.  and  the  non-linearity 

set  at  6  V,  the  dead-beat  controller  is 

1  -  0.697Z1  -  0.178Z2  +  0.1046Z3 

D (Z)  =  4.087  - — r — 

1  -  0.05Z1  -  0.675Z2  -  0.275Z 

This  system  was  then  simulated  on  the  analog  computer  to  study 
the  effect  of  changes  of  parameters  on  the  output.  The  system  was 


i 
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observed  to  behave  very  well.  The  various  disturbances  studied  are 
shown  in  table  3.1. 

Step  input  to  the  system  =  2  V 

Settling  time  refers  to  the  time  required  for  the  output  to  return 
to  its  normal  value  after  the  application  of  the  disturbance. 

The  small  blips  on  the  curves  show  the  sampling  instants. 

TABLE  3.1 


Fig. 

Parameter 

Original 

New 

Settling 

Changed 

Value 

Value 

Time 

3.2 

- 

1 

Response  without 

disturbance 

3.2 

- 

2 

k2 

1 

1.1 

2T 

3.2 

- 

3 

k2 

1 

1.5 

5T 

3.2 

- 

4 

k2 

1 

2.0 

7T 

3.2 

_ 

5 

Time  constant 

1 

1 

4T 

S(S  +  1) 

S(S  +  1.1) 

3.2 

_ 

6 

Time  constant 

1 

1 

5T 

S(S  +  1) 

S(S  +  1.5) 

3.2 

7 

Time  constant 

1 

1 

7T 

S(S  +  1) 

S(S  +  2) 

3.2 

- 

8 

A  step  disturbance  of  5  V  is 

applied  at  the 

output 

14T 

STEP  INPUT 

TO  THE  SYSTEM 

=  5  V 

3.3 

- 

1 

Response  without 

disturbance 

3.3 

- 

2 

k2 

1 

2.0 

8T 

3.3 

3 

Time  constant 

1 

1 

13T 

S(S  +  1) 

S(S  +  2) 

3.3 

- 

4 

A  step  disturbance  of  5  V  is 

applied  at  the 

output  11T 


■ 


Original 


Fig.  Parameter 

Changed 


New 

Value  Value 


Settling 

Time 


STEP  INPUT  TO  THE  SYSTEM  =  8  V 


3.4 

-  1 

Response  without  disturbance 

3.4 

-  2 

K2  1 

2.0 

8T 

3.4 

-  3 

Time  constant  1 

1 

9T 

S(S  +  1) 

S(S  +  2) 

3.4 

-  4 

A  step  disturbance  of  5  V  is 

applied  at  the 

output 

12T 

STEP  INPUT  TO  THE  SYSTEM 

=  10  V 

3.5 

-  1 

Response  without  disturbance 

3.5 

-  2 

K2  1 

2.0 

9T 

3.5 

-  3 

K2  1 

0.5 

6T 

3.5 

-  4 

Time  constant  1 

1 

4T 

S(S  +  1) 

S(S  +  l.|5) 

3.5 

-  5 

Time  constant  1 

1 

4T 

S(S  +  1) 

S(S  +  1.2) 

3.5 

-  6 

A  step  disturbance  of  5  V  is 

applied  at  the 

output 

19T 

When 

the  plant  was  changed  to  1 

,  the  disturbanc 

S(S  + 

1.25) 

could  not  be  removed.  The  error  kept  on  increasing  in  the  branch 
of  the  controller^ 

It  may  be  appropriate  now  to  make  the  following  comments  about 
the  controller. 

1.  The  controller  was  designed  for  ramp  disturbances  only.  It  was 
discussed  in  chapter  two  that  a  change  in  parameter  is  equival¬ 
ent  to  the  sum  of  exponential,  step  and  ramp  disturbances  at 


■ 
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Fig.  3.3 
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Fig.  3.5 
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the  output. 


The  controller  was  designed  on  the  assumption  that  there  was  no 
other  input  to  the  non-linearity  except  that  due  to  the  ramp 
disturbance.  When  the  parameter  was  changed,  there  was  already 
a  signal  at  the  input  to  the  non-linearity.  This  signal  was 
due  to  the  step  input  to  the  system. 


The  controller  was  designed  for  the  plant 


,  but 


S(S  +  1) 

when  the  parameter  was  changed  the  controller  had  to  deal  with 
a  different  plant. 

The  results  in  the  table  3.1  show  that  the  controller  gives 
a  satisfactory  performance  under  the  changed  conditions. 

Next  the  step  disturbance  was  applied  at  different  locations 
Fig.  3.6  :  Step  input  to  the  system  =  2  V 

3.6  -  1  :  Response  without  disturbance 

3.6  -  2  :  A  step  disturbance  of  5  V  is  applied  at  the  point 

'1'  in  fig.  3.1 

3.6  -  3  :  A  step  disturbance  of  5  V  is  applied  at  the  point 

'm'  in  fig.  3.1. 

The  saturation  limit  of  the  non-linearity  was  then  varied.  The 
results  are  shown  in  table  3.2. 

TABLE  3.2  (a) 

The  saturation  limit  of  the  plant  non-linearity  is  changed  to 


5  V. 
Fig. 


Parameter 

Changed 


Original 

Value 


New 

Value 


Settling 

Time 


3.7  -  1 


Response  without  disturbance 


. 


> 

' 


. 
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Fig. 

Parameter 

Original 

New 

Settling 

Changed 

Value 

Value 

Time 

3.7  -  2 

A  step  disturbance  of  5  V 

is  applied  at 

the 

output 

13T 

3.7  -  3 

k2 

1 

2.0 

8T 

3.7  -  4 

Time  constant 

1 

1 

8T 

S(S  +  1) 

S(S  +  2) 

3.7  -  5 

k2 

1 

0.5 

6T 

3.7  -  6 

Time  constant 

1 

1 

6T 

S(S  +  1) 

S(S  +  0.5) 

TABLE  3.2  (b) 

The 

saturation  limit 

of  the  plant  : 

non-linearity  is  changed  to 

7  V. 

THE  STEP 

INPUT  TO  THE  SYSTEM  =  2  V 

Fig. 

Parameter 

Original 

New 

Settling 

Changed 

Value 

Value 

Time 

3.8  -  1 

Response  without 

disturbance 

3.8  -  2 

A  step  disturbance  of  5  V 

is  applied  at 

the 

output 

18T 

3.8  -  3 

k2 

1 

0.5 

7T 

3.8  -  4 

k2 

1 

2.0 

8T 

3.8  -  5 

Time  constant 

1 

1 

9T 

S(S  +  1) 

S(S  +  2) 

3.8  -  6 

Time  constant 

1 

1 

S(S  +  1) 

S(S  +  0.5) 

When  the  plant  was  changed  to  _ 1 _  the  oscillations 

S(S  +  0.5) 

observed  at  the  output  exhibited  a  slow  decay. 


Table  3.1  shows  that  the  settling  time  for  the  disturbance 
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Fig.  3.6 


Fig.  3.7 
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Fig.  3.8 


Fig.  3.10 
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increases  as  the  step  input  to  the  system  is  increased.  This  is 
because  of  the  fact  that  with  an  increase  in  the  step  input  to  the 
system,  the  output  of  the  non-linearity  approaches  the  saturation 
limit.  The  non-linearity  is  then  unable  to  give  large  corrective 
outputs  required  to  remove  disturbances. 

When  the  input  to  the  system  was  10  V  (signal  at  the  non¬ 
linearity  was  5  V) ,  the  controller  could  not  remove  a  change  of 
more  than  207,  in  the  time  constant.  The  corrective  signal  re¬ 
quired  was  larger  than  the  margin  available  between  the  operating 
point  and  the  saturation  limit  (1  V  in  this  case) .  The  error  kept 
on  increasing  at  the  input  to  the  controller. 

A  decrease  in  time  constant,  say  to  1  ,  was  still 

S(S  +  0.5) 

tolerated.  The  corrective  signal  required  in  this  case  was 
negative.  This  was,  however,  achieved  at  the  cost  of  a  more  oscil¬ 
latory  response.  The  (controller  produces  both  positive  and  negative 
steps  during  the  transient  response^ whereas  the  negative  steps 
were  reproduced  at  the  output  of  the  non-linearity,  the  positive 
steps  were  limited  in  magnitude  to  the  margin  available  between  the 
operating  point  and  the  saturation  limit.  Thus  if  the  system  is 
operating  near  the  saturation  limit,  a  disturbance  requiring  a 
larger  output  from  the  non-linearity  may  not  be  removed  while  a 
disturbance  requiring  a  smaller  output  will  be  removed  with 
oscillations  and  in  a  longer  settling  time. 

If  a  ramp  input  is  applied  to  the  above  system?the  non-linearity 
will  soon  reach  its  saturation  limit  and  then  the  above  discussion 
applies.  The  system  is  thus  not  very  good  for  ramp  inputs  or  for 
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those  step  inputs  which  make  the  non-linearity  to  operate  near  its 
saturation  limit.  This  will  be  true  for  a  conventional  feedback 
system  also  and  is  thus  not  a  drawback  of  the  conditional  feedback 
system  alone. 

3.2  SATURATION  TYPE  NON-LINEARITY  IN  BETWEEN  TWO  TRANSFER  FUNCTIONS 


The  saturation  type  non-linearity  was  then  placed  as  shown  in 
fig.  3.9. 


K2  =  1 

Fig.  3.9  T  =  1  sec. 


The  dead-beat  controller  was  designed  by  the  rep-op  method.  In 
a  case  like  this,  the  design  of  a  dead-beat  controller  by  the  state 
variable  approach  will  be  quite  involved  but  it  was  quite  simple 
by  the  rep-op  method. 

The  dead-beat,  controller  for  a  ramp  disturbance  of  2  V  and  the 
saturation  limit  of  the  non-linearity  set  at  6  V  was  found  to  be 

1  -  1.01.3Z1  +  0.252Z2 

D(Z)X  =  3.922  - - — 

1  -  0,5  Z1  -  0.5Z2 

T  =  1  sec. 

This  controller  was  then  modified  by  trial  to  give  a  good  re- 
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sponse  to  a  step  input  also.  The  coefficients  of  Z1  and  Z2  in 

the  numerator  were  varied.  The  various  controllers  studied  were  as 

follows 

1  -  1.1Z1  +  0.260Z2 

D(Z)2  =  3.922  - - r - - - - 

1  -  0.5Z1  -  0.5Z2 

1  -  1.15Z1  +  0.270Z2 
D(Z)3  =  3.922  — - - 

1  -  0.5Z1  -  0.5Z2 

1  -  1.2Z1  +  0.280Z2 

D(Z)4  =  3.922  - - 

1  -  0.5Z1  -  0.5Z2 

1  -  1.125Z1  +  0.265Z2 

D(Z)  5  =  3.922  - - - - - 

1  -  0.5Z1  -  0.5Z2 

The  results  are  shown  in  table  3.3. 


TABLE  3.3 


Controller 

Max.  Overshoot 

to  a  step 

input 

Settling  Time 

for  a  step 

input 

Settling  Time 

for  a  ramp 

input 

D(Z)  x 

80% 

10T 

3T 

D(Z)2 

60% 

11T 

9T 

D(Z)  3 

42% 

14T 

12T 

D(Z)4 

30% 

28T 

20T 

D(Z)  5 

50% 

13T 

10T 

D(Z)5  was  finally  chosen. 

This  system  was  then  simulated  on  the  analog  computer  to  study 
the  effect  of  change  of  parameters  on  the  output.  The  results  are 


shown  in  table  3.4 


. 
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TABLE  3.4 


Fig. 

Parameter 

Original 

New 

Settling 

Changed 

Value 

Value 

Time 

STEP  INPUT 

TO  THE  SYSTEM  =  2  V 

3.10 

- 

1 

Response  without  disturbance 

3.10 

- 

2 

k2 

1 

1.1 

- 

3.10 

- 

3 

k2 

1 

1.5 

6T 

3.10 

- 

4 

k2 

1 

2.0 

9T 

3.10 

- 

5 

k2 

1 

0.7 

6T 

3. 10 

- 

6 

k2 

1 

0.5 

7T 

3.11 

- 

1 

Response  without  disturbance 

3.11 

- 

2 

Time  constant 

1 

1 

- 

S(S  + 

1)  S(S 

+  1.1) 

3.11 

— 

3 

Time  constant 

1 

1 

8T 

S(S  + 

1)  S  (S 

+  1.5) 

3.11 

- 

4 

Time  constant 

1 

1 

8T 

S(S  + 

1)  S(s 

+  2) 

3.11 

— 

5 

Time  constant 

1 

1 

7T 

S(S  + 

1)  S(S 

+  0.7) 

3.11 

• 

6 

Time  constant 

1 

1 

12T 

S(S  + 

1)  S(S 

+  0.5) 

STEP  INPUT 

TO  THE  SYSTEM  =  10  V 

3.12 

- 

1 

Response  without  disturbance 

3.12 

- 

2 

k2 

1 

1,1 

- 

3.12 

- 

3 

k2 

1 

2.0 

11T 

3.12 

- 

4 

k2 

1 

0.5 

6T 

3.12 

5 

Time  constant 

1 

1 

3T 

S(S  +  1)  S(S  +  1.1) 
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Fig. 

Parameter 

Original 

New 

Settling 

Changed 

Value 

Value 

Time 

3.12  -  6 

Time  constant 

1 

1 

9T 

S(S  +  1) 

S(S  +  2) 

3.12  -  7 

Time  constant 

1 

1 

14T 

S(S  +  1) 

S(S  +  0.5) 

STEP  INPUT 

TO  THE  SYSTEM  = 

16  V 

3.13  -  1 

Response  without 

disturbance 

3.13  -  2 

k2 

1 

1.2 

7T 

3.13  -  3 

k2 

1 

2.0 

13T 

3.13  -  4 

k2 

1 

0.5 

8T 

3.13  -  5 

Time  constant 

1 

1 

5T 

S(S  +  1) 

S(S  +  1.2) 

3.13  -  6 

Time  constant 

1 

1 

9T 

S(S  +  1) 

S(S  +  2) 

3.13  -  7 

Time  constant 

1 

1 

11T 

S(S  +  1) 

S(S  +  0.5) 

It  was  observed  in  this  case  that  it  was  possible  to  apply  step 
inputs  which  make  the  non-linearity  operate  in  its  saturation 
region,  but  only  the  parameter  variations  before  the  non-linearity 
will  be  tolerated.  The  output  of  the  model  is  6  V/sec,  This  is 
also  the  maximum  value  of  the  output  required  from  the  plant. 
Therefore  as  long  as  the  plant  non-linearity  has  a  saturation  limit 
of  6  V  or  slightly  more,  the  plant  will  always  be  able  to  yield  an 
output  of  6  V/sec.  In  fact  those  parameter  variations  which  tend  to 
produce  larger  signals  at  the  input  to  the  non-linearity  may  not 
be  able  to  affect  the  output  at  all.  The  output  of  the  non¬ 
linearity  will  stay  at  6  V.  This  is,  of  course,  true  for  a 
perfect  non- linear ity#A  s  will  be  seen  later,  if  the  saturation 
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Fig.  3.12 
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Fig.  3.13 


Fig.  3.14 
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limit  of  the  plant  non-linearity  is  less  than  that  of  the  model 
non-linearity,  then  it  would  not  be  possible  to  apply  step  inputs 
which  saturate  the  plant  non-linearity. 

The  response  to  step  disturbances  is  shown  in  table  3.5. 


TABLE  3.5 


STEP  INPUT  TO  THE  SYSTEM  =  2  V 


Fig. 
3.14  -  1 
3.14  -  2 

3.14  -  3 

3.14  -  4 


Remarks 


Response  without  disturbance 
A  step  disturbance  of  5  V  is  applied  at 


the  output 

A  step  disturbance  of  5  V  is  applied  at 
the  point  1 1'  in  fig.  3.9 
A  step  disturbance  of  5  V  is  applied  at 
the  point  'm'  in  fig.  3.9.  A  disturbance 
at  'm'  gives  rise  to  an  equal  disturbance 
at  the  output 

Table  3.6  shows  the  response  when  the  saturation  limit  of  the 
plant  non-linearity  was  changed  to  5  V. 


TABLE  3.6 

STEP  INPUT  TO  THE  SYSTEM  -  2  V 


Fig. 

Parameter 

Original 

New 

Settling 

Changed 

Value 

Value 

Time 

3.15  -  1 

Response 

without 

disturbance 

3.15  -  2 

k2 

1 

1.2 

- 

3.15  -  3 

k2 

1 

0.7 

4T 

3.15  -  4 

Time  constant 

1 

1 

8T 

S(S  +  1)  S(S  +  2) 
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Fig. 

Parameter 

Original 

New 

Settling 

Changed 

Value 

Value 

Time 

3.15  - 

5 

Time  constant 

1 

1 

3T 

S(S  +  1) 

S(S  +  1.2) 

3.15  - 

6 

Time  constant 

1 

1 

6T 

S(S  +  1) 

S(S  +  0.7) 

3.15  - 

7 

A  step  disturbance  of  5  V  is 

applied  at 

the  output 

11T 

A 

step 

input  of  8  V  was 

then  applied  to 

the  system  and 

the 

above  parameter  variations  were  studied.  The  system  behaved  well 
in  that  case  except  that  with  large  parameter  changes,  the 
settling  time  was  larger. 

A  step  input  of  16  V  will  not  work  in  this  case,  because  the 
model  will  be  giving  a  ramp  output  of  6  V/sec.  and  since  the  satura¬ 
tion  limit  in  the  plant  non-linearity  is  5  V,  the  plant  output  can¬ 
not  exceed  5  V/sec.  The  error  will  therefore  increase  continu¬ 
ously  in  the  branch  of  the  controller. 

The  saturation  limit  of  the  non-linearity  was  then  changed  to 
7  V.  The  responses  to  parameter  variations  are  shown  in  table  3.7. 

TABLE  3.7 

STEP  INPUT  TO  THE  SYSTEM  =  2  V 


Fig. 

Parameter 

Original 

New 

Settli 

Changed 

Value 

Value 

Time 

3.16  -  1 

Response  without 

dis  turbance 

3.16  -  2 

k2 

1 

1.2 

- 

3.16  -  3 

k2 

1 

1.5 

7T 

3.16  -  4 

k2 

1 

0.7 

5T 
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Fig. 

Parameter 

Original 

New 

Settling 

Changed 

Value 

Value 

Time 

3.16  -  5 

Time  constant 

1 

1 

4T 

S(S  +  1) 

S(S  +  1.2) 

3.16  -  6 

Time  constant 

1 

1 

13T 

S(S  +  1) 

S(S  +  0.7) 

3.16  -  7 

A  step  disturbance  of  5 

V  is  applied  at 

the  output 

10T 

STEP  INPUT  TO 

THE  SYSTEM  = 

20  V 

3.17  -  1 

Response  without  disturbance 

3.17  -  2 

k2 

1 

1.2 

6T 

3.17  -  3 

k2 

1 

2.0 

16T 

3.17  -  4 

k2 

1 

0.5 

8T 

3.17  -  5 

Time  constant 

1 

1 

3T 

S(S  +  1) 

S(S  +  1.2) 

3.17  -  6 

Time  constant 

1 

1 

8T 

S(S  +  1) 

S(S  +  2) 

3.17  -  7 

Time  constant 

1 

1 

14T 

S(S  +  1) 

S(S  +  0.5) 

3.17  -  8 

A  step  disturbance  of  20 

V  is  applied  at 

the  output 


In  this  case  when  the  step  input  is  applied  to  the  system, 
the  output  of  the  plant  tends  to  approach  7  V/sec.  and  that  of 
the  model  to  6  V/sec.  The  controller  starts  acting  and  brings 
the  output  of  the  plant  to  6  V/sec.  This  means  that  the  plant  is 
still  operating  in  the  linear  range  of  the  non-linearity  and  there 
is  some  margin  available  for  additional  signals  to  correct  the 
disturbances.  The  settling  time  will,  however,  be  larger  than  for 
the  case  when  a  smaller  step  input  is  applied  to  the  system. 


? 
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Fig.  3.16 
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Fig.  3  o 17 
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3.3  HYSTERESIS  TYPE  NON-LINEARITY 

This  type  of  non-linearity  is  often  encountered  in  the 
sensing  devices.  The  hysteresis  non-linearity  was,  therefore, 
placed  in  the  feedback  path.  The  model  was  simulated  similar  to 
the  plant  i.e.  it  contained  a  sampler,  and  a  zero  order  hold, 
but  no  non-linearity  was  included  in  it.  It  was  found  that  the 
model  without  a  non-linearity  worked  as  good  as  the  model  with  a 
non-linearity.  It  was,  therefore,  decided  not  to  include  the  non¬ 
linearity  in  the  model.  This  will  be  true  in  practice  because  of 
economic  reasons. 

The  dead-beat  controller  used  here  was  the  same  as  in  article 

2.1  for  the  plant  1  .  The  conditional  feedback  arrangement  is 

S  +  1 


shown  in  fig.  3.19. 


Fig.  3.19 


T  =  1  sec, 


The  above  system  was  then  simulated  on  the  analog  computer  to 
study  the  response  to  a  step  input  and  also  the  sensitivity  of  the 
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system  to  parameter  variations  and  external  disturbances.  The 
simulation  of  hysteresis  non-linearity  is  described  in  appendix  B. 

The  effect  of  hysteresis  was,  that  the  output  deviated  from  its 
normal  value  by  half  the  value  of  the  hysteresis  as  shown  in  fig. 3. 20, 


Fig.  3.20 


On  the  computer,  the  response  seemed  to  drift  between  its 
upper  and  lower  limits.  The  limits  are  shown  in  fig.  3.20  (a). 
This  is  due  to  the  noise  produced  in  the  computer,  but  in  actual 
practice  the  output  would  stay  at  one  limit  or  the  other. 

The  various  disturbances  studied  are  shown  in  table  3.8. 


TABLE  3.8 

STEP  INPUT  TO  THE  SYSTEM  =  10  V 


Fig. 

Parameter 

Changed 

Original 

Value 

3.18  -  1&2 

k2 

1 

3.18  -  3 

k2 

1 

3.18  -  4&5 

k2 

1 

3.21  -  1 

k2 

1 

K2 


New 

Value 

1.02 

1.01 

1.005 

1.1 


3.21  -  2 


1 


1.3 
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Fig. 


3.21  -  3 
3.21  -  4 

3.21  -  5 

3.22  -  1 
3.22  -  2 
3.22  -  3 
3.22  -  4 

3.22  -  5 

3.23  -  1 
3.23  -  2 
3.23  -  3 

3.23  -  4 

3.24  -  1 

3.24  -  2 

3.25  -  1 
3.25  -  2 
3.25  -  3 

3.25  -  4 

3.26  -  1 
3.26  -  2 
3.26  -  3 
3.26  -  4 

3.26  -  5 


Parameter 

Changed 

K2 

k2 

k2 

k2 

k2 

k2 

k2 

k2 

Ki 

Kl 

Ki 

Ki 

Ki 

Ki 

Ki 

Ki 

Ki 

Ki 

Ki 

Ki 

Ki 

Time  constant 


Original 

Value 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

(S  +  1) 

1 

(S  +  1) 


New 

Value 

1.5 

0.95 

0.99 

0.98 

0.95 

0.90 

0.70 

0.50 

1.005 

1.01 

1.02 

1.05 

1.1 

1.5 

0.995 

0.99 

0.98 

0.95 

0.90 

0.70 

0.50 

1 

(S  +  1.01) 
1 

(S  +  1.05) 


Time  constant 


-  63 


Fig. 

Parameter 

Original 

New 

Changed 

Value 

Value 

3.27  -  1 

Time  constant 

1 

1 

(S  +  1) 

(S  4-  1.11) 

3.27  -  2 

Time  constant 

1 

1 

(S  +  1) 

(S  +  1.2) 

3.27  -  3 

Time  constant 

1 

1 

(S  +  1) 

(S  +  1.3) 

3.27  -  4 

Time  constant 

1 

1 

(S  +  1) 

(S  +  1.5) 

3.27  -  5 

Time  constant 

1 

1 

(S  +  1) 

(S  +  2) 

3.28  -  1 

Time  constant 

1 

1 

(S  +  1) 

(S  +  0.99) 

3.28  -  2 

Time  constant 

1 

1 

(S  +  1) 

(S  +  0.95) 

3.28  -  3 

Time  constant 

1 

1 

(S  +  1) 

(S  +  0.9) 

3.28  -  4 

Time  constant 

1 

1 

(S  +  1) 

(S  +  0.7) 

When  the  plant  was  changed 

to  1 

,  the  output 

(S  +  0.5) 

became  very  oscillatory.  This  is  due  to  the  fact  that  the 


controller  had  been  designed  for  a  plant  1  .  If  there 

(S  +  1) 

was  a  large  veriation  in  the  plant,  the  controller  did  not  act 
effectively. 

The  above  observations  show  that  any  disturbance  which  is  able 
to  produce  a  signal  at  the  output  of  the  hysteresis  will  soon  be 
detected  by  the  controller  and  its  effect  on  the  output  removed. 

A  ramp  inputvas  then  applied  to  the  system.  The  response  was 
good  as  long  as  there  was  no  variation  in  parameter.  A  parameter 
variation  would  give  rise  to  a  ramp  disturbance  at  the  output  in 
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Fig.  3.24 
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Fig.  3.25 


6  7" 


Fig 
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the  steady  state.  The  controller  was  designed  for  step  disturbances 
only.  This  controller  would  give  a  constant  error  to  ramp  dis¬ 
turbances.  The  constant  error  will  appear  at  the  input  to  the 
controller.  Since  the  controller  contains  a  free  digital  in¬ 
tegrator,  its  output  will  keep  on  increasing.  This  is  exactly 
what  was  observed  on  the  computer. 

This  drawback  can,  however,  be  eliminated  if  the  controller  is 
designed  for  both  ramp  and  step  inputs. 

A  continuous  data  model  could  also  be  used  here  because  no 
non-linearity  is  associated  with  the  model.  This  was  tried.  The 
continuous  data  model  used  was  the  same  as  derived  in  article  2.1. 
The  response  was  as  good  as  with  the  sampled  data  model.  The 
continuous  data  model  will  be  economical  to  use  in  practice.  But 
the  use  of  a  continuous  data  model  will  be  possible  when  the  non¬ 
linearity  is  present  only  after  the  sensing  device  and  the  forward 
path  transfer  functions  are  linear. 
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CHAPTER  IV 


The  conditional  feedback  can  also  be  used  to  optimise  the  response 

of  oscillatory  systems.  The  dead-beat  controller  is  designed  as  usual 

depending  upon  the  types  of  disturbances  expected.  The  model  is  then 

modified  so  that  the  response  of  the  oscillatory  system  is  improved. 

A  second  order  oscillatory  system  has  a  transfer  function  of 

the  form  _ U)^1 _  ?  <1  |  *  The  conditional  feedback 

i  UHv2- 

arrangement  for  this  is  shown  in  fig.  4.1. 


Fig.  4.1 

The  above  system  was  simulated  on  the  analog  computer.  It  was  found 
that  the  response  of  the  oscillatory  system  could  be  improved  by 
setting  the  damping  ratio  of  the  model  to  a  different  value  from  that 
of  the  plant.  The  optimum  response  for  a  plant  with  the  transfer 

v 

function  _ 1 _  was  found  by  trial.  This  is  shown  in 

S2  +  1.2S  +  1 

table  4. 1 

5(p  =  damping  ratio  of  the  plant 
=  damping  ratio  of  the  model 

Settling  time  refers  to  the  time  required  for  the  response  to  decrease 
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to  and  stay  within  21  of  its  final  value. 


STEP  INPUT  TO  THE  SYSTEM  =  10  V 
Sampling  period  T  =  0.5  sec. 

Dead-beat  controller  found  by  rep-op  method 


"l 

1  -  0.1234Z 

D  (Z)  =  1.931  - — 

1  -  0.425Z 


"2 

1.243Z  +  0.7335Z 


0.4Z  -  0.175Z 


0.6 


TABLE  4.1 


Fig. 

5  'VYN 

Max.  overshoot 

or  undershoot 

Settling  time 

4.2  - 

1 

1 

1% 

9  sec . 

4.2  - 

2 

0.8 

1 . 5% 

7  sec . 

4.2  - 

3 

0.7 

2 . 51 

3.6  sec. 

4.2  - 

4 

0.75 

21 

4.4  sec . 

4.2  - 

5 

0.72 

11 

3.6  sec. 

4.3  - 

1 

0.73 

11 

3 . 8  sec . 

4.3  - 

2&3 

0.74 

21 

4  sec. 

4.3  - 

4 

0.75 

21 

4.4  sec . 

4.3  - 

5 

0.  72 

1% 

3.6  sec . 

The  above  data  indicates  an  optimum  value  of  ?  to  be  0.72. 

J'm 

The  effect  of  sampling  period  of  the  controller  on  the  optimum 
response  was  then  studied.  It  was  found  that  as  the  sampling  period 
is  reduced,  the  response  improves.  The  disadvantage  of  reducing  the 
time  period  is  that  the  controller  is  required  to  give  larger  steps. 
If  these  large  steps  saturate  the  plant,  it  will  not  be  possible  to 


, 
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reduce  the  sampling  period  beyond  a  certain  point.  The  effects  of 


variation  of  sampling  time  period  on  the  response  are  shown  in  table 
4.2.  All  dead-beat  controllers  were  designed  by  the  rep-op  method. 


1.  Transfer  function  of  plant  =  _ 1 _ 

S2  +  S  +  1 


h  -  °-5 

T  =  1  sec. 


"1  "2 

1  -  0.861Z  +  0.394Z 

D (Z)  =  1.852  - - 

"1  "2 
1  -  0.6Z  -  0.4Z 


T  =  0. 5  sec. 


"1  "2 

1  -  0.1043Z  -  1.293Z  +  0.812Z 

D (Z)  =  2.333  - - r— 

“l  "2 

1  -  0.234Z  -  0.516Z  -  0.25Z 


T  =  0. 2  sec. 


"1  “2 

1  -  0.7730Z  -  1.115Z  4-  0.899Z 

D  (Z)  =  13.94  — - ^ - 

“1  “2 

1  -  0.25Z  -  0.5Z  -  0.25Z 


2.  Transfer  function  of  plant  =  1 

S2  +  0.6S  +  1 


%  -  0.3 

T  =  1  sec. 


"1  “2 

1  -  0.9175Z  +  0.5775Z 

D  (Z)  =  1.523  - r - - 

1  -  0.5625Z1  -  0.4375Z2 


T  =  0 . 5  sec . 
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"1  "2 

1  -  0.1124Z  -  1.484Z  +  1.079Z 

D(Z)  =  2.035  - r - 

“1  “2 

1  -  0.25Z  -  0.5Z  -  0.25Z 


T  =  0.2  sec . 


"i  “2  "3 

1  -  0.7723Z  -  1.175Z  -  Z 

D (Z)  =  13.23  - . 

"1  "2  ~ J 

1  -  0.25Z  -  0.5Z  -  0.25Z 


TABLE  4.2 


Fig. 

4 

T  Max. 

Settling 

( 

overshoot 

time 

JYr\ 

*4.4 

-  1 

0.5 

14% 

8  sec. 

- 

4.4 

-  2 

0.5 

1  sec.  47o 

7.8  sec . 

0.663 

4.4 

-  3 

0.5 

0.5  sec.  27> 

6.4  sec . 

0.667 

4.4 

-  4 

0.5 

0.2  sec.  17o 

3.8  sec. 

0.73 

*4.4 

-  5 

0.3 

33% 

12  sec. 

- 

4.4 

-  6 

0.3 

1  sec.  167. 

11.2  sec. 

0.555 

4.4 

-  7 

0.3 

0.5  sec.  127. 

10.8  sec. 

0.55 

4.4 

-  8 

0.3 

0.2  sec.  3% 

8  sec . 

0.6375 

The 

effect  of 

change  of  parameters  on 

the  output  were 

then 

*  These  curves  show  the  response  of  the  transfer  function 


_  (jo>x _  to  a  step  input,  i.e.  the  model  and 

(Sd-2)  (  S*>  2-i^Sd-  toy*) 

the  controller  are  out  of  circuit.  The  other  curves  show  the 
improvement  of  this  response  by  the  use  of  conditional 


feedback 
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Fig.  4,3 
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Fig.  4.4 


Fig.  4.5 
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studied.  The  results  are  shown  in  table  4.3. 

STEP  INPUT  TO  THE  SYSTEM  =  50  V 

The  plant  is  _ 1 _ 

S2  +  S  +  1 

The  dead-beat  controller  is 

"1  "2 

1  -  0. 7730Z  -  1. 115Z  +  0.899Z 

D(Z)  =  13.94  - - - - 

“1  “2 

1  -  0.25Z  -  0.5Z  -  0.25Z 

.Iyyy  *  0 .  73 .  T ~  0>  Z.  ^ec. 


The  parameters  refer  to  fig.  4.1 

TABLE  4.3 


Fig. 

Parameter 

Original 

New 

Changed 

Value 

Value 

4.5-1 

Kl 

1 

1.1 

4.5  -  2 

Ki 

1 

1.3 

4.5  -  3 

Ki 

1 

1.5 

4.5  -  4 

Ki 

1 

0.9 

4.5  -  5 

Kl 

1 

0.7 

4.5  -  6 

Kl 

1 

0.5 

4.6  -  1 

k2 

1 

1.  1 

4.6  -  2 

k2 

1 

1.3 

4.6  -  3 

k2 

1 

1.5 

4.6  -  4 

k2 

1 

0.9 

4.6  -  5 

k2 

1 

0.7 

4.6  -  6 

k2 

1 

0.5 

The  effects 

of  change  of  the 

natural  frequency  and 

the  damp 

ratio  of  the  plant  on  the  transient  response  are  shown  in  tables 


.  V 
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4. 4  and  4.5, 

The  optimum  response  from  the  oscillatory  plant  can  also  be 
achieved  by  a  dead-beat  controller  and  the  conventional  feedback. 

The  conventional  feedback  will,  in  fact,  be  economical  to  use  because 
it  does  not  require  a  model.  The  conditional  feedback  was,  there¬ 
fore,  compared  with  conventional  feedback  and  it  was  found  that  the 
former  is  superior  to  the  latter  in  performance.  The  conditional 
feedback  system  was  much  less  sensitive  to  parameter  changes  then 
the  conventional  feedback  system.  In  addition  to  this  the  conditional 
feedback  system  has  other  advantages  which  were  discussed  in 
chapter  I. 

The  two  systems  were  compared  on  the  following  two  bases: 

1.  Same  sampling  period  for  the  controllers. 

2.  Same  rise  time  of  the  response  to  a  step  input. 

The  comparison  was  made  on  the  basis  of  transient  response. 

In  the  first  case  (same  sampling  period  for  the  controllers) , 
the  response  to  a  disturbance  will  be  similar  for  both  types  of  feed¬ 
back,  because  the  feedback  loop  for  disturbances  is  essentially  the 
same.  The  response  to  disturbances  will  be  better  for  conditional 
feedback  in  the  second  case  (same  rise  time  to  step  inputs) .  The 
disturbances  will  be  removed  much  quicker  and  with  less  overshoot, 
because  of  the  shorter  sampling  period.  The  arrangement  of  the  two 
types  of  systems  is  shown  in  fig.  4.7. 

1 .  The  same  sampling  period  for  the  controllers 

Sampling  period  -  0.2  sec. 


77 


"1  "2 

1  -  0.7988Z  -  0.9270Z  +  0.798Z 

D (Z)  =  14.25  - o - 

“1  ”2  “ J 

1  -  0.25Z  -  0.5Z  -  0.25Z 


This  controller  was  used  for  both  the  types  of  feedback  systems. 


STEP  INPUT  TO  THE  SYSTEM  =  50  V 


(a)  Conditional  feedback 


(b)  Conventional  feedback 


Fig.  4.7 


Conventional  feedback  Conditional  feedback 
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Fig.  4.11 


Fig.  A. 12 
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2 .  The  same  rise  time  of  the  response  to  a  step  input 

Sampling  period  of  controller  for  conditional  feedback  =0.2  sec. 


D,(Z) 

"1  "2 

(1  -  0.7988Z  -  0.9270Z  4-  0.798Z  ) 

=  14.25  - t - 

"1  "2 

(1  -  0.25Z  -  0.5Z  -  0.25Z  ) 

Rise  time  of  response  to  a  step  input  =  4  sec. 

Sampling  period  of  controller  for  conventional  feedback  =  1.33  sec 


D  (Z) 

“1  “2 

1  +  0. 113Z  +  0. 1395Z 

=  0.813  - - 

“I  "2 

1  -  0.662Z  -  0.338Z 

Rise  time 

of  response  to  a  step  input  =  4  sec. 

The  results  are  shown  in  table  4.5 


Step  input  to  the  system  =  50  V 
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Fig.  4.13 


Fig.  4. 14 


\ 


-  87 


Fig.  4.15 

These  results  clearly  show  that  the  conditional  feedback 
system  is  much  less  sensitive  to  parameter  changes  than  the  con¬ 
ventional  feedback  system. 
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CONCLUS ION 

This  thesis  indicates  that  the  idea  of  conditional  feedback 
can  be  advantageously  applied  to  linear  and  non-linear  sampled  data 
systems.  The  feedback  loop  can  be  designed  so  that  the  disturbances 
have  an  appreciably  small  effect  on  the  output.  The  input-output 
response  is  independent  of  this  feedback  loop.  The  nature  of 
the  input-output  response  can  be  entirely  different  from  the  dis¬ 
turbance-  output  response. 

Non-linearities  can  be  tolerated  in  the  system  without  causing 
instability  with  reference  to  the  input  signal.  The  response  of 
oscillatory  systems  can  also  be  optimised  by  manipulating  the  model. 
The  inclusion  of  a  dead-beat  controller  into  the  conditional  feed¬ 
back  configuration  does  not  make  the  output  highly  sensitive  to 
parameter  changes  as  it  would  do  in  the  conventional  feedback 
configuration. 
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APPENDIX  A 

REPETITIVE  OPERATION  METHOD  FOR  DESIGNING  THE  DEAD-BEAT  CONTROLLERS 

(REP -OP  METHOD) 


This  is  an  analog  computer  method  -for  designing  the  dead-beat  con¬ 
trollers  for  linear  and  non-linear  systems.  This  method  was  developed 
by  Professor  Y.  J.  Kingma. 

Fig.  A. 1  shows  a  closed  loop  control  system. 


D(Z) 

I -  - 1 


The  controller  D(S)  is  designed  to  give  a  dead-beat  response.  When 
the  input  is  applied  to  the  system,  the  digital  controller  produces  a 
series  of  step  inputs  such  that  the  response  of  the  system  is  optimum. 
The  duration  of  each  step  input  is  T.  The  number  of  these  steps  de¬ 
pends  upon  the  order  of  the  system,  the  type  of  the  system  and  the  type 
of  the  input  to  the  system.  The  digital  controller  can  be  expressed 
in  the  z-domain  as  follows. 

"l  “2 

ho  +  hj_Z  4-  h2Z  4-  h3Z  4- - 

D(Z)  =  - - - =i - =3 - 

eo  4-  e^Z  4-  e2Z  4-  e3Z  4- - 

If  this  is  analyzed  by  the  state  Transition  Technique,  it  is  found 


that  ho,  hl5  h2,  I13 


etc.  are  the  magnitudes  of  the  step  inputs 


—  <  -f 
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to  the  plant  at  t  =  0,  t  =  T,  t  =  2T,  t  =  3T  -  etc.  sampling 

instants  respectively  and  eO,  el9  e2,  63  -  etc.  are  the  mag¬ 

nitudes  of  the  signals  that  are  fed  into  the  controller  at  the  sam¬ 
pling  instants  t  =  0,  t  =  T,  t  =  2T,  t  =  3T  -  etc.  respectively. 

Thus  if  the  magnitudes  of  the  step  inputs  into  the  plant  and  the  mag¬ 
nitudes  of  the  signals  into  the  controller  at  the  sampling  instants 
are  known  the  dead-beat  controller  can  be  written  down.  This  infor¬ 
mation  can  be  obtained  by  the  rep-op  method.  The  procedure  is 
explained  below. 

The  transfer  function  of  the  plant  is  simulated  on  the  analog 
computer  without  the  feedback  loop.  The  computer  is  set  in  the  re¬ 
petitive  operation  mode.  The  input  to  the  system  is  applied  from  the 
special  equipment  available  for  this  purpose.  This  equipment  is  shown 
in  fig.  A. 2.  The  box  shown  in  fig.  A. 2  gives  a  series  of  steps  as 
shown  in  fig.  A. 3.  The  magnitude  and  the  sign  of  each  of  these  steps 
is  controlled  by  a  separate  potentiometer  and  a  switch.  The  output 
of  the  plant  transfer  function  at  various  locations  is  displayed  on  a 
multi-channel  oscilloscope.  The  outputs  appear  stationary  on  the 
screen  because  the  computer  is  in  the  high  speed  repetitive  operation 
mode.  The  potentiometers  and  the  switches  are  then  adjusted  till  all 
the  outputs  meet  the  requirements  of  dead-beat  response.  The  settings 
of  the  potentiometers  and  the  sign  of  the  switches  give  the  value  of 
the  h's,  and  the  values  of  the  e's  can  be  read  off  on  the  oscilloscope. 
The  dead-beat  controller  can  then  be  written  down. 

The  box  in  fig.  A. 2  can  give  five  steps,  the  maximum  value  of 
each  step  is  11.0  volts.  The  sampling  period  can  be  adjusted  by  the 


t  e  * c  t 
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Fig .  A0 2 
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potentiometer  marked  T  in  fig.  A. 2.  The  magnitudes  of  the  steps  is 
adjusted  by  the  potentiometers  numbered  1  to  5. 

This  method  is  simple,  convenient  and  direct.  There  is  no  inter¬ 
sample  ripple  in  the  response.  The  presence  of  a  non-linearity  or 
a  transportation  delay  does  not  make  the  method  any  more  difficult. 

This  method  is  also  very  convenient  to  design  dead-beat  controllers  for 
plants  which  have  auxiliary  feedback  loops.  The  method  is  explained 
for  the  following  cases; 

A. 1  A  third  order,  type  one  system  with  a  step  input 

Plant  transfer  function  =  1  . 

S(S  +  1)  (S  +  2) 

The  computer  diagram  is  shown  in  fig.  A. 4. 


0.2 

Fig.  A. 4 


The  minimum  number  of  sampling  periods  required  for  dead¬ 
beat  response  in  this  case  is  three.  The  form  of  the  signals 
at  channels  A,  B,  C  and  D  is  shown  in  fig.  A. 5.  At  t  =  3T, 
the  signal  at  channel  A  should  reach  the  value  of  the  step 
input  and  after  that  it  should  stay  at  that  value.  The  signals 
at  channels  B  ancf  C  should  reduce  to  zero  at  t  -  3T  and  stay  at 


zero  after  that. 
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(c)  Channel  C 


Fig.  A. 5 

The  values  found  in  this  case  were 

+50 
-25. 1 
+2.5 

The  dead-beat  controller  will  be 

“1  -4 

1  -  0.502Z  +  0.05Z 


eo  = 

14.41 

ho  = 

el  = 

9.88 

hi  = 

e2  = 

1.15 

h2  = 

D  (Z)  =  3.47 


1  +  .685Z1  +  0.0798Z^ 
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A. 2  A  second  order,  type  one  system  containing  a  saturation  type 

non-linearity  in  between  the  transfer  functions.  A  ramp  input 

is  applied  to  the  system 

The  computer  diagram  is  shown  in  fig.  A. 6.  It  is  assumed 

here  that  the  ramp  input  will  be  applied  to  the  system  at  one  of 

the  sampling  instants. 

Plant  transfer  function  =  1 

S(S  +  1) 


The  signals  at  the  various  channels  are  shown  in  fig.  A. 7. 

No  signal  will  appear  at  the  input  to  the  plant  at  the 
first  sampling  instant  because  the  input  is  a  ramp  and  there  is 
a  sampler  and  a  zero  order  hold  before  the  plant. 

In  this  case  it  is  required  that  the  controller  should  also 
act  as  an  integrator  i.e.  it  should  hold  the  value  at  the  output 
when  the  input  to  it  reduces  to  zero.  This  is  necessary  because 
the  integrator  in  the  plant  will  give  a  ramp  output  only  when  it 
has  a  constant  step  input  to  it.  When  the  ramp  output  equals 
the  ramp  input  (dead-beat  response) ,  the  signal  in  the  error 
branch  will  be  zero  and  so  will  be  the  input  to  the  controller. 
This  integrating  effect  can  be  achieved  by  associating  a  digital 


% 
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Channel  A 


Channel  B 


Channel  C 


Fig.  A. 7 
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integrator  of  the  form  shown  in  fig.  A. 8  with  the  digital  con¬ 
troller  . 


The  digital  integrator  can  be  expressed  as  1  This  will 

(i  -  z1) 

be  multiplied  into  the  D(Z)  found  by  the  above  method. 

The  requirements  of  dead-beat  response  in  this  case  are  as 
follows 

Channel  A;  -  The  signal  at  channel  A  should  acquire  the  value  of 
the  ramp  input  at  t  =  3T  and  then  stay  at  that  value  after  that. 
Channel  B:  -  The  signal  at  channel  B  should  acquire  a  constant 
value  at  t  =  3T  and  stay  at  that  value  after  that.  If  the  ramp 
input  is  k  V/sec.  then  the  constant  value  is  k. 

Channel  C:  -  In  this  case  the  last  step  applied  to  the  system 
should  be  of  magnitude  k,  if  the  ramp  input  is  k  V/sec.  This 
last  step,  in  the  above  case,  will  be  applied  at  t  =  3T.  This 
is  so,  because  at  t  =  3T,  the  input  to  the  controller  will  reduce 
to  zero  and  then  the  digital  integrator  will  hold  the  value  of 
the  last  step.  This  last  step  will  be  an  additional  step.  The 
error  will  have  already  been  reduced  to  zero  at  t  =  3T.  The 


- 
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information  obtained  will  be  of  the  following  form 


In  a  digital  controller,  the  powers  of  the  numerator  and 
the  denominator  have  to  be  equal.  Thus  in  the  above  case  e2 
is  required.  This  e2  will  be  obtained  when  the  digital  integrator 
is  multiplied  into  the  digital,  controller  obtained  by  the  above 
insufficient  information.  The  modified  controller  will  be 

“1  "2 
hQ  +  h±Z  +  h2Z 

D  (2)  m  = 

(eQ  +  exZ  )(1  -  Z~)  A.  1 

The  digital  integrator  inserted  between  the  dead-beat  con¬ 
troller  and  the  plant  will  modify  the  steps  being  produced  by  the 
controller.  This  modification  can  be  explained  as  follows. 

From  equation  A. 1,  the  steps  produced  by  the  digital  controller 
are  hQ,  hx  and  h2.  The  digital  integrater  is  shown  in  fig.  A. 8. 
When  the  step  inputs  to  this  integrator  are  hG,  h2  and  h2  at 
the  instants  t.  =  0,  t  =  T  and  t  =  2T  respectively,  the  outputs 
of  the  integrator  will  be 


ho 

at 

t  =  0 

A. 2 

ho 

+ 

hl 

at 

t  =  T 

A. 3 

^o 

+ 

hi  4-  h2 

at 

t  =  2T 

A.  4 

But  to  get  the  dead-beat  response  from  the  plant,  the  step  inputs 
into  the  plant  should  be  hQ,  h1  and  h,2  at  the  instants  t  =  0, 
t  -  T  and  t  -  2T  respectively.  This  can  be  achieved  by  modifying 
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the  step  outputs  of  the  controller.  Let  the  modified  step 
outputs  of  the  controller  be 
at  t  =  0  -+  hQ ' 
at  t  =  T  — >  h-L ' 
at  t  =  2T  — > h2 1 

Solving  equations  A. 2,  A. 3  and  A. 4  for  hD  1 ,  h-;  '  and  h2 1  yields, 

^o  '  =  ^o 

hi'  =  hl  "  ho 
h2 '  =  h2  -  hx 

Therefore  the  modified  dead-beat  digital  controller  will  be 


"1  “2 
ho'  +  hx'Z  +  h2'Z 

D(Z)m  =  ~i  “i 

(eo  +  &iZ  )  ( 1  -  Z  ) 


h0'  +  hx'Z  +  h2 ' Z 

“1  "2 

eO  +  (el  _eo)  ^  +  el^ 


It  may  be  noted  here  that  the  sum  of  the  coefficients  of  the 
powers  of  Z  in  the  denominator  will  always  be  equal  to  zero. 

For  the  plant  considered  above  the  values  found  were* 
ramp  input  =  20  V/sec. 

saturation  limit  of  the  non-linearity  =  60  V 


eo  =  +20 


ex  =  +10 


e2  =  0 


Hence 


hQ  =  +78.43 


h,  =  -1.03 


h2  =  +20 


h0 '  =  +78.43  h2'  =  +  21.03 

h-L '  =  -79.46 


-  100 


“1  "2 

78.43  -  79.46Z  +  21.03Z 

D(Z)m  =  — - — -  ■ 

(20  +  10Z  ) (1  -  Z  ) 


“1  "2 

78.43  -  79.46Z  +  21.03Z 


“1  "2 
20  -  10Z  -  10Z 


A. 3  A  second  order  type  zero  oscillatory  system,  with  a  step  input 

The  plant  transfer  function  1 

S2  +  S  +  1 

The  requirements  in  this  case  are  exactly  the  same  as  for 
the  type  one  system  with  a  ramp  input.  The  controller  should  act 
as  an  integrator  in  this  case  also.  The  computer  diagram  is  shown 
in  fig.  A. 9. 


Fig.  A. 9 


The  requirements  for  dead-beat  response  are 

1.  The  signal  at  channel  A  should  reach  the  value  of  the  step 
input  at  t  =  2T  and  remain  at  that  value  after  that. 

2.  The  signals  at  B  and  C  should  reduce  to  zero  at  t  =  2T 
and  then  remain  at  zero  after  that. 

3.  The  step  output  of  the  controller  at  t  =  2T  should  be  equal 


-  ; 
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to  the  step  input  to  the  system. 

The  values  found  were 

eQ  =  3.0  hQ  =  +86,16 

ei  =  1.47  h-L  =  -66.54 

e2  =  0  h2  =  3.0 

Hence 

hG'  =  +86. 18 
h±'  =  -152.72 
h2 '  =  +69.54 


D(Z)m 


“1  "2 

86.18  -  152. 72Z  +  69.54Z 

"l  ~i 

(3  +  1.47Z  ) (1  -  Z  ) 


86.18 


3  - 


“1  "2 
-  152. 72Z  +  69.54Z 


"1  "2 
1.53Z  -  1.47Z 


Thus  it  may  be  seen  that  it  is  very  convenient  and  straight 
forward  to  design  the  digital  controllers  by  the  rep-op 
method. 

A  digital  computer  program  was  worked  out  to  check 
the  correctness  of  the  dead-beat  controllers  designed  by 
the  rep-op  method.  This  program  is  shown  in  appendix  C. 
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APPENDIX  B 

This  appendix  illustrates  the  simulation  of  the  hysteresis  type 
non-linearity  on  the  analog  computer.  This  scheme  minimizes  the 
drift.  The  arrangement  is  shown  in  fig.  B.l. 


Fig.  B.l 


APPENDIX  C 


Digital  computer  program  to  check  the  correctness  of  dead-beat 

controllers 

This  digital  computer  program  gives  the  values  of  the  output 
at  the  sampling  instants  if  the  coefficients  of  the  powers  of  Z 
in  the  Z-transforms  of  the  plant  and  the  digital  controller  are 
known. 

A  step  input  to  the  system  is  considered. 

The  closed  loop  system  is  shown  in  fig.  C.l. 


Fig.  C.l 


The  Z-transforms  of  the  controller  and  the  plant  can  be  written  as 

D4Z4  +  D3Z3  +  D2Z2  +  DXZ  +  D0 

D(Z)  =  P-l  - - 

4  3 

E4Z  +  E3Z  +  E2Z2  +  EXZ  +  E0 

4  3^ 

G4Z  +  G3Z  +  G2Z2  +  GXZ  +  Go 

Gp(Z)  =  P2  - 

^  4  3 

H4Z  +  H3Z  +  H2Z2  +  HXZ  +  H0 

D(Z)  Gp(z) 


C(Z) 


1  +  D(Z)  Gp(z) 


x  R(Z) 


V 
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R(Z)  =  Z 

Z  -  1 

10  9  8 

A(I)  Z  +  A(I  -  1)Z  +  A(I  -  2)Z  + . 

/.  C(Z)  =  - - - - 

10  9  8 

B  ( I)  Z  +  B  ( I  -  1)Z  +  B  ( I  -  2)  Z  + . 

1.  The  values  of  D'S,  E'S,  G'S,  H'S,  P'S  are  to  be  read  into  the 

program  after  the  end  card.  The  value  of  I  is  also  to  be  read. 

If, 

m  =  Max.  power  of  Z  in  the  denominator  of  D(Z) 

n  =  Max.  power  of  Z  in  the  denominator  of  Gp(Z) 

then  I  =  m  +  n  +  2  (Maximum  value  of  I  for  this  program  is  10) 

The  program  is  as  under 

DIMENSION  Z(100),  B(20)  ,  A(20) 


READ 

(5,1) 

DO, 

Dl, 

D2, 

D3, 

D4 

READ 

(5,1) 

E0, 

El, 

E2, 

E3 , 

E4 

READ 

(5,1) 

GO, 

Gl, 

G2, 

G3, 

G4 

READ 

(5,1) 

HO, 

HI, 

H2 , 

H3, 

H4 

READ 

(5,2) 

I 

READ 

(5,3) 

PI, 

P2 

1  FORMAT  (IX,  5 (F10 . 7 ,  IX)) 

2  FORMAT  (IX,  12) 

3  FORMAT  (IX,  2(F10.6,  IX)) 

A(  1)  =  +  0.0000000 

A(2)  =  P1*P2*D0*G0 

A(3)  =  P1*P2*(D1*G0  +  D0*G1) 

A (4)  =  P1*P2*(D2*G0  +  D1*G1  +  D0*G2) 

A(5)  =  Pl*P2*(D3*G0  +  D2*G1  +  D1*G2  +  D0*G3) 
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A(6)  =  P1*P2*(D4*G0  +  D3*G1  +  D2*G2  +  D1*G3  +  D0*G4) 
A(  7)  =  P1*P2*(D4*G1  +  D3*G2  +  D2*G3  4-  D1*G4) 

A(8)  =  P1*P2*(D4*G2  +  D3*G3  +  D2*G4) 

A(9)  =  P1*P2*(D4*G3  +  D3*G4) 

A(  10)  =  P1*P2*D4*D4 
CO  =  EO*HO  +  A(2) 

Cl  =  (El*HO  +  E0*H1  +  A(3) ) 

C2  =  E2*H0  +  E1*H1  +  E0*H2  +  A(4» 


C3  =  (E3*H0 
C4  =  (E4*H0 
C5  =  (E4*H1 
C6  =  (E4*H2 
C7  =  (E4*H3 


+  E2*H1  +  E1*H2  +  E0*H3  +  A(5)) 

+  E3*H1  +  E2*H2  +  E1*H3  +  E0*H4  +  A(6)) 
+  E3*H2  +  E2*H3  +  El*H4  +  A(7)) 

+  E3*H3  +  E2*H4  +  A(8)) 

+  E3*H4  +  A(9) ) 


C8  =  (E4*H4  +  A(  10) ) 


B  ( 1)  =  -CO 

B (2)  =  CO  -  Cl 

B (3)  =  Cl  -  C2 

B (4)  =  C2  -  C3 

B(5)  =  C3  -  C4 

B (6)  =  C4  -  C5 

B ( 7)  =  C5  -  C6 

B (8)  =  C6  -  C7 

B (9)  =  C7  -  C8 

B  ( 10)  =  C8 
ZO  =  A(I)/B(I) 

Z ( 1)  =  ( A ( I  -  1)  -  ZO*B(I  -  1))/B(I) 
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40  J  =  I  -  1 

DO  60  K  =  2, J 
Z  (KK)  =0.0 
JJ  =  K  -  1 
DO  61  N  =  1,JJ 
MM  =  I  -  K  +  N 

61  Z  (KK)  =  Z  (KK)  +  Z(N)*b(MM) 

LLL  =  I  -  K 

60  Z  (K)  =  (A(LLL)  -  ZO*B(LLL)  -  Z(KK))/B(I) 
50  L  =  I  -  1 

DO  20  K  =  1,35 

LL  =  L  +  K 
Z  (LL)  =0.0 
II  =  I  -  1 
DO  21  N  =  1,11 
M  =  N  -  1  +  K 

21  Z  (LL)  =  Z  (LL)  +  Z(M)*d(N) 

20  Z  (LL)  =  -  (Z)  LL)  )  /B  ( I) 

WRITE  (6,4)  ZO 
DO  25  I  =  1,43 
25  WRITE  (6,4)  Z(I) 

4  FORMAT  (5X, F10 . 5) 

WRITE  (6,6) 

6  FORMAT  ( 1H1 , 5X,  14HVALUES  OF  A  (I)) 

DO  26  I  =  1 , N 


26  WRITE  (6,4)  A(I) 


; 
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WRITE  (6,7) 

7  FORMAT  ( 1H1 , 5X, 14HVALUES  OF  B(I)) 

DO  27  I  =  1,10 
27  WRITE  (6,4)  B(I) 

5  FORMAT  ( 10X,  F10 . 5) 

END 

Example 

A  second  order  oscillatory  system  with  the  transfer 

function  1 _  was  chosen 

S2  +  S  +  1 

Z  +  0.9251 

Gp(Z)  =-0.0187  — - - - 

Z2  -  1.783Z  4-  0.819 

The  dead-beat  controller  found  by  the  rep-op  method  was 

“l  "2  "3 

1  -  0.7988Z  "  0.9270Z  +  0.798Z 

D(Z)  =  14.25  - - - - - - - 

-I  -2  -3 

1  -  0.25Z  -  0.5Z  -  0.25Z 

T  -  0 . 2  sec . 

1  =  7 

Input  to  the  system  =  a  unit  step  U(t) 

The  results  found  were  as  follows: 


P.T.O. 


VALUES  OF  THE  OUTPUT  AT  THE  SAMPLING  INSTANTS 


t  =  0 
t  =  T 
t  =  2T 
t  =  3T 
t  =  4T 
t  *  5T 
t  -  6T 
t  =  7T 
t  =  8T 
t  =  9T 
t  «  10T 
t  =  11T 
t  =  12T 
t  =  13T 
t  =  14T 
t  =  15T 
t  =  16T 
t  «  17T 
t  =  18T 
t  =  19T 

A(  1) 

A(2) 

A(3) 

A  (4) 
A(5) 


0.00000 

0.26647 

0.77087 

1.00519 

1.00238 

1.00520 

1.00603 

1.00546 

1.00326 

1.00208 

1.00103 

0.99940 

0.99883 

0.99840 

0.99811 

0.99795 

0.99789 

0.99792 

0.99802 

0.99818 

0.00000 

0. 19665 

-0.01595 

-0.44394 

0.03366 


VALUES  OF 


t  =  20T 
t  =  2  IT 
t  =  22T 
t  =  23T 
t  =  24T 
t  =  25T 
t  =  26T 
t  =  27T 
t  =  28T 
t  a  29T 
t  =  30T 
t  =  3  IT 
t  =  32T 
t  =  33T 
t  =  34T 
t  =  35T 
t  a  36T 
t  a  37T 
t  a  38T 
t  =  39T 

MU 

A(6) 

A(7) 

A(8) 

A(9) 

A(  10) 


0.99837 

0.99859 

0.99882 

0.99905 

0.99927 

0.99948 

0.99966 

0.99983 

1.00009 

1.00018 

1.00025 

1.00029 

1.00032 

1.00033 

1.00033 

1.00031 

1.00029 

1.00026 

1.00023 

1.00019 

0.26647 

0.00000 

0.00000 

0.00000 

0.00000 
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VALUES  OF  B  ( I) 


B  ( 1) 

0.00810 

B  (6) 

-2.76652 

B  (2) 

-0.02840 

B  (  7) 

1.00000 

B  (3) 

0.02749 

B  (8) 

0.00000 

B  (4) 

-0.80560 

B(9) 

-0.00000 

B  (5) 

2.56493 

B  ( 10) 

0.00000 

